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ABSTRACT 


This  investigation  was  undertaken  with  the  objective  of  developing 
techniques  for  solving  the  problem  of  sound  transmission  from  a 
harmonic  monopole  source  through  a  finite,  corrugated  boundary  between 
fluid  media  and  to  corroborate  these  techniques  by  means  of  laboratory 
data  obtained  for  the  important  case  of  air-to-water  sound  transmission* 

Using  a  form  of  the  Kirchoff  radiation  integral,  expressions  were 
derived  describing  both  tie  refected  and  transmitted  acousHc 
potentials  due  to  the  scattering  of  monopcle  sound  at  a  corrugated 
boundary  between  fluids.  For  the  case  of  sound  transmission  through 
the  boundary  the  integral  expression  was  reduced  to  an  approximate 
closed-form  solution  which  is  valid  in  the  limit  of  gecmetncal  optics 
and  for  extended  interface  geometries  for  which  the  tangent  plane 
approximation  is  applicable. 

Empirical  data  was  obtained  using  four  different  laboratory  models 
of  corrugated  (sinusoidal)  surface  sections  designed  to  float  on  the 
surface  of  an  anechoic  tank  so  as  to  simulate  different  corrugated  air- 
water  boundaries.  Numerous  comparisons  of  theoretical  results  with 
experimental  data  obtained  by  loudspeaker  insonif ication  of  these 
model  surfaces  demonstrated  the  suitability  of  the  theory  and  illu¬ 
strated  the  inaccuracies  which  arise  in  cases  for  whLch  the  tangent 
plane  approximation  is  unsuitable. 
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CHAPTER  I 


INTRODUCTION 

General  Introduction  to  the  Problem  Area 

In  recent  years,  a  considerable  amount  of  interest  has  evolved 
regarding  changes  which  occur  in  the  characteristics  of  acoustic 
signals  which  are  transmitted  within  or  between  two  media  of  different 
characteristic  impedances  which  share  a  common  boundary.  Of 
particular  interest  is  the  problem  of  predicting  the  influence  of 
various  rough  boundary  surface  configurations  on  the  characteristics 
of  these  signals.  Until  very  recently,  such  studies  have  been  almost 
entirely  devoted  to  the  problem  of  sound  scattering  (reflection)  from 
such  boundaries  by  various  classes  of  rough  surfaces.  Evidently, 
this  interest  in  acoustic  signal  scattering  was  motivated  principally 
by  the  need  to  understand  the  effects  of  rough  surfaces  such  as  that 
of  the  sea  surface  or  ocean  bottom  on  various  echo-ranging  signals 
transmitted  in  the  ocean.  However,  with  the  development  of  more 
sophisticated  underwater  acoustics  equipment,  it  has  become  obvious 
that  reception  of  acoustic  signals  passing  between  air  and  water  is 
practicable.  Hence,  in  order  to  more  fully  exploit  this  capability, 
renewed  interest  has  been  shown  in  studies  related  to  the  problem 
of  sound  transmission  through  various  classes  of  rough  surfaces 
separating  fluid  media.  Moreover,  in  addition  to  problems  involving 
sound  transmission  between  different  fluids  separated  by  rough 
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boundaries,  various  related  problems  abound  involving  the  transmission 
of  vibrational  waves  between  fluid-solid  and  solid-solid  interfaces 
in  such  fields  as  oceanography  and  seismology.  Hence,  basic  studies 
which  enhance  our  knowledge  of  vibrational  wave  transmission  between 
media  separated  by  arbitrary  boundary  geometries  can  be  expected  to 
contribute  significantly  to  the  advancement  of  technology  even  in 
diciplines  not  directly  involved  with  the  transmission  of  sound  waves, 
per  se. 

Need  for  the  Study 

Present  theoretical  literature  contains  but  scant  reference  to 
the  problem  of  sound  transmission  through  a  rough  interface  separating 
two  media  of  different  characteristic  impedances.  That  which  does 
exist  deals  in  various  approximate  ways  with  the  transmission  of  a 
perfectly  planar  incident  wave  at  the  surface  of  an  infinite  boundary 
described  as  a  general  periodic  uneven  surface  or  a  boundary*  comprised 
of  a  corrugated  surface  extending  to  infinity  in  two  directions.  No 
reference  whatever  has  been  found  which  relates  to  the  point-source 
rough  interface  transmission  problem.  Evidently,  one  of  the  reasons 
that  this  latter  problem  has  been  virtually  neglected  to  date  is  that 
the  mathematical  difficulties  encountered  in  the  development  of  a 
rigorous  theoretical  treatment  have  proven  to  be  intractable.  In  fact, 
Lysanov  (1958)  points  out  that  it  is  Impossible  to  use  the  method  of 
separation  of  variables  to  obtain  a  solution  of  the  wave  equation 
which  satisfies  the  boundary  conditions  specified  on  the  uneven 
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surface.  Denied  the  use  of  this  powerful  method  for  solving  the 
relevant  system  of  equations,  it  is  not  difficult  to  appreciate  the 
mathematical  complexities  which  result. 

For  the  specific  case  of  a  simple  (monopole)  source  in  one  medium, 
both  (wave)  integral  equations  and  geometrical  op tics -type  (ray  theory) 
solutions  exist  which  can  be  used  to  compute  the  sound  field  developed 
in  the  adjacent  medium,  but  only  for  the  case  of  a  smooth  (flat) 
boundary  geometry.  Much  atii?  remains  to  be  done  to  extend  these 
techniques,  or  develop  new  ones,  which  will  permit  one  to  compute  the 
sound  field  in  the  adjacent  media  for  rough  boundary  geometries. 

Finally,  no  empirical  data  of  any  consequence  appears  to  be 
available  for  the  case  of  sound  transmitted  between  fluid  media 
separated  by  rough  boundaries. 

Specific  Statement  of  the  Problem 

The  objective  of  this  study  is  to  advance  our  understanding  of 
the  processes  involved  in  acoustic  transmission  between  fluid  media 
which  are  separated  by  an  uneven  boundary  surface.  In  this  context, 
the  specific  problem  herein  chosen  for  study  concerns  the  description 
of  the  steady-state  sound  field  which  arises  in  a  fluid  medium  due  to 
the  transmission  of  acoustic  energy  from  a  point  (monopole)  harmonic 
source  located  in  a  second  (adjacent)  fluid  medium  sharing  a  common, 
limited  time- invariant  corrugated  boundary  with  the  first. 
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Scope  and  Limitations  of  the  Study 

In  order  to  advance  our  understanding  of  acoustic  transmission 
between  media  separated  by  rough  boundary  surfaces  this  study  was 
initiated  to: 

(1)  Develop  a  mathematical  model  and  computational  procedure  to 
predict  the  acoustic  field  in  a  fluid  medium  due  to  the  operation  of 
a  point  (monopole)  harmonic  source  in  a  second  fluid  medium  of 
different  acoustic  impedance.  Each  of  the  two  homogeneous  isotropic 
media  are  assumed  to  occupy  its  own  infinite  half-space.  Transmission 
of  acoustic  energy  is  assumed  to  take  place  only  through  a  finite 
portion  of  the  common  boundary  between  the  two  media.  Further,  the 
cross-section  of  the  acoustically  transparent  interface  between  the 
media  is  assumed  to  be  composed  of  regular  corrugations  (sinusoidal 
undulations) . 

(2)  Obtain  empirical  data  on  the  transmission  of  sound  from  air 
to  water  through  a  finite  corrugated  boundary  surface  with  excitation 
frequency  and  source  an.;le  as  parameters. 

(3)  Establish  the  suitability  of  the  mathematical  model  by 
comparing  predicted  results  with  the  empirical  data  obtained  during 
the  laboratory  tests. 

(4)  Investigate  various  relationships  between  the  physical 
parameters  of  the  surface  and  source  with  the  objective  of  defining 
the  conditions  under  which  the  accuracy  of  the  theoretical  procedures 


will  be  maintained. 
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To  facilitate  the  manipulation  and  evaluation  of  the  expressions 
to  be  developed  herein,  it  was  necessary  to  introduce  certain  specific 
assumptions  in  the  course  of  the  mathematical  analysis.  Actually, 
these  restrictions  have  become  familiar  in  studies  of  this  type  and 
most,  if  not  all  of  them,  have  been  used,  either  implicitly  or 
explicity,  in  all  of  the  studies  on  scattering  and  refraction  reported 
to  date.  These  limitations  include  the  following  suppositions: 

1.  Multiple  reflections  in  either  medium  can  be  ignored. 

2.  The  source  and  receiver  are  restricted  to  positions  not 

in  the  immediate  vicinity  of  the  interface  in  the  sense 
that  kjRg  and  are  large. 

3.  Masking  (shadowing)  of  one  portion  of  the  surface  by 
adjacent  ones  need  not  considered. 

4.  Transient  effects  are  not  included. 

3.  The  tangent  plane  approximation  is  valid. 

6.  The  spatial  rate  of  change  of  the  phase  of  the  incident 

wave  changes  much  more  rapidly  than  that  of  the  transmission 
or  reflection  coefficients  in  the  region  cf  the  interface. 

Note  that  no  specific  restriction  on  the  height  of  the  surface 
is  made,  providing  the  above  assumptions  remain  valid. 

Previous  Related  Studies 

Apparently,  the  earliest  and  certainly  the  best  known  approximate 
approach  to  solving  the  problem  of  reflection  and  transmission  of 
harmonic  waves  at  a  rough  interface  is  that  described  by  Rayleigh  in 


his  "Theory  of  Sound"  published  in  1877  (See  Reference  8).  For  a 
normally  incident  plane  harmonic  wave  and  an  infinite  sinusoidal 
surface,  Rayleigh  formulated  a  solution  to  the  problem  in  terms  of 
an  equivalent  problem  in  which  the  solution  of  an  infinite  set  of 
linear  equations  must  be  found.  Inverting  this  system,  he  obtained 
a  solution  to  the  transmitted  wave  in  the  form  of  an  infinite  series 
made  up  of  individual  terms  representing  distinct  planar  waves  (the 
spectra)  which  are  scattered  in  various  directions.  In  order  to 
evaluate  the  coefficients  in  the  series,  it  is  necessary  to  assume 
that  the  acoustic  wavelength  is  much  smaller  than  the  wavelength  of 
the  surface  undulations,  observer  and  source  are  in  the  far  field  of 
the  surface  and  the  slope  of  the  uneven  surface  is  small. 

A  noteworthy  extension  of  Rayleigh’s  basic  treatment  has  been 
reported  by  Asano  (1966).  He  assumes  that  a  general  periodic  surface 
can  be  decomposed  into  a  Fourier  series  composed  of  elementary 
corrugations  and  then  proceeds  to  fit  boundary  conditions  at  a 
representative  corrugated  interface,  much  in  the  same  way  as  Rayleigh. 
However,  in  addition  to  considering  purely  longitudinal  vibrations, 
he  includes  transverse  vibrations  (shear  waves)  in  his  formulation  for 
several  cases  of  interest  in  seismology.  Although  his  work  is,  for 
the  most  part,  restricted  to  the  case  of  plane  normally  incident 
elastic  waves,  he  does  include  the  subject  of  oblique  incidence  when 
the  incident  wave  is  of  the  plane  longitudinal  (P)  type.  Again,  to 
facilitate  application  of  the  boundary  conditions,  both  the  amplitude 
and  slope  of  the  corrugations  were  assumed  to  be  small. 
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By  numerically  evaluating  the  first  two  coefficients  of  the 
resulting  series  solution,  Asano  reached  several  important  conclusions 
regarding  the  effects  of  corrugation  amplitude  and  wavelength  on  the 
strengths  of  the  lower  order  spectral  components.  The  most  striking 
observation  noted  was  the  presence  of  large,  rapid  fluctuations  in 
the  values  of  the  series  coefficients  as  the  ratio  of  the  corrugation 
wavelength  to  the  wavelength  of  the  incident  P-wave  ACA.)  was  varied 
in  the  range  of  about  one-half  to  two.  Especially  rapid  changes  were 
frequently  encountered  for  values  of  AgA^)  just  slightly  less  than 
one.  On  the  other  hand,  for  values  of  AgAj)  in  excess  of  two,  the 
amplitude  coefficients  tended  to  remain  reasonably  constant.  Although 
his  numerical  evaluation  was  limited  to  velocity  contrasts  less  than 
two  to  one,  Asano  noted  a  tendancy  toward  greater  fluctuation  in  the 
amplitudes  of  all  low  order  coefficients  as  the  velocity  contrast  was 
increased.  He  also  noted  a  tendency  for  the  regularly  refracted 
P-wave  to  gain  energy  at  the  expense  of  the  regularly  reflected  P-wave 
as  the  corrugation  amplitude  was  increased.  For  the  case  of  an 
obliquely  directed  incident  P-wave,  the  angle  of  incidence  appeared  to 
have  little  effect  on  the  regularly  reflected  and  refracted  P-waves 
but  a  somewhat  more  significant  (although  still  moderate)  influence 
on  the  irregularly  reflected  and  refracted  P-waves. 

A  particularly  illuminating  study  of  a  related  problem  in  the 
field  of  electricity  and  magnetism  is  that  of  Park  and  Erteza  (1969). 
Undertaking  to  solve  the  problem  of  radiation  from  an  arbitrarily 
oriented  electromagnetic  dipole  source  in  the  presence  of  a  rough, 
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Imperfectly  conducting  earth,  they  present  a  form  of  the  vector 
Helmholtz  integral  which  is  suited  for  their  application.  By  means  of 
this  formulation,  they  proceed  to  express  the  reflected  and  transmitted 
Hertz  potentials  in  terms  of  the  dyadic  reflection  and  transmission 
coefficients  of  the  boundary.  Applying  themselves  to  a  mathematical 
model  of  a  slightly  rough  surface  described  by  a  random  process  with 
a  stationary  Gaussian  height  distribution,  they  obtain  solutions  of 
the  geometrical  optics  type  for  the  expected  Hertz  potentials  and 
expected  power  arising  from  steady  state  excitation  of  an  elementary 
dipole.  In  addition  to  being  a  worthwhile  exposition  of  an  application 
of  the  mathematics  of  geometrical  optics  to  transmission  in  the 
presence  of  a  rough  boundary,  this  recent  contribution  is  apparently 
the  first  significant  work  to  be  reported  on  the  problem  Gf  scattered 
electromagnetic  transmission  through  a  rough  interface. 

Mashashvili  and  Urusovskii  (1969)  report  an  extension  of  work 
previously  reported  by  Urusovskii  (1965)  in  which  they  consider  the 
problem  of  plane  wave  diffraction  by  an  extended  periodically  uneven 
interface  between  two  homogeneous  media.  Using  the  Green rs  formula, 
together  with  the  Hankel  function  representation  of  the  incident  wave 
field,  they  make  use  of  the  periodicity  of  the  integral  kernels  to 
develop  a  solution  for  the  reflected  and  transmitted  pressure  fields 
in  terms  of  a  set  of  algebraic  equations. 


CHAPTER  I 


THEORETICAL  DEVELOPMENT 
Kirchoff  Radiation  Integral 

Due  to  the  similarity  in  the  development  of  the  integral 
equations  for  scattering  (reflection)  and  transmission  (refraction) 
arising  from  point  source  insonification  of  a  corrugated  boundary,  it 

is  useful  to  carry  through  the  general  theoretical  development  for 

* 

both  cases.  However,  methods  for  evaluating  the  resulting  integral 
expressions  as  well  as  application  of  corroborative  experimental  work 
will  be  restricted  to  the  case  of  transmission  through  the  boundary,  as 
indicated  previously.  However,  it  is  important  to  emphasize  that  only 
minor  modifications  are  required  to  adapt  the  computational  procedures 
developed  for  sound  transmission  to  the  case  of  sound  scattering 
(reflection) . 

An  integral  equation  for  the  portion  of  the  incident  sound  field 
scattered  back  into  the  original  medium  by  the  corrugations,  as  well 
as  an  expression  for  the  field  transmitted  into  the  adjacent  medium, 
may  be  found  with  the  aid  of  the  Kirchoff  integral.  In  order  to 
develop  a  form  of  the  Kirchoff  integral  suited  for  use  in  the  present 
study,  we  begin  by  consideration  of  Green’s  Theorem  in  the  form 
(Reference  4,  Page  15): 

fv  =  jfg  (*  H  -  y  |*)ds,  (2.1) 
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where  +  and  i|r  are  arbitrary  scalar  fields,  dv  is  the  differential 
element  of  volume,  and  ds  the  differential  element  of  surface  area. 
The  partial  derivative  denotes  differentiation  with  respect  to  the 
unit  normal  directed  outward  from  the  volume  under  consideration. 

For  the  solution  of  the  Hyperbolic  (wave)  equation  relevant  to 
the  present  case,  we  require  specification  of  Cauchy  boundary 
conditions  (the  potential  and  its  normal  derivative)  on  an  "open” 
surface  (Reference  4,  Page  17).  Here,  the  open  surface  is  the  three- 
dimensional  half-space  at  a  fixed- time.  We  can  manipulate  Equation 
(2.1)  to  obtain  an  integral  expression  for  the  velocity  potential  at 
any  point  within  the  particular  medium  to  be  used  in  Equation  (2.1) 
and  for  any  time,  t,  greater  than  the  fixed  initial  instant  chosen. 
However,  it  is  expedient  to  restrict  ourselves  to  the  steady-state 
time  periodic  case  from  the  outset  and,  with  this  understanding,  we 
can  develop  the  pertinent  form  of  the  Kirchoff  equation  with  minimal 
effort.  Refer  to  the  geometry  depicted  by  Figure  1.  Then  using  the 
subscript,  i,  to  denote  the  numerals  1  or  2  (which  refer  to  medium  1 
or  2)  and  the  subscript,  j,  to  indicate  the  letters  r  or  t  for 
reflection  or  transmission,  we  let 


*  = 


e-i(cot-kiR.) 


R 


i 


(2.2) 


Separating  the  time  variation  from  the  potential,  we  can  write 


j 


-io)t 

e 


(2.3) 


Now,  che  physics  of  Che  problem  requires  ChaC  Che  scaCCered 


poCenCials  satisfy  a  homogeneous  wave  equation,  i.e., 

or,  using  Equation  (2.3), 

v2*j +  ki2  = 0  :  ki2  £  72 


(2.4) 


(2.5) 


Further,  since  in  Equation  (2.2)  is  Green’s  function  for  the  medium 
containing  the  observation  point,  it  satisfies  the  equation 


^  -  h  ~2t  =  ‘4n  B(Ri)  e'1<ut 

c  ^  dt 


(2.6) 


or 


2  e^i^i  2  e^^i^i  — 

> +  ki  (7  >  =  *4*  S<V- 


ik.R,  „  aik,R 

*  K 

Using  Equations  (2.2),  (2.3),  (2.5)  and  (2.7),  it  is  easy  to  show 
that  Equation  (2.1)  becomes: 

ik  R, 


(2.7) 


„  i  n  Ow #  i  _ 

*J  <V  V  zi>  =  kfs  Tt  +  iki  <!  +  ^j]ds-  (2-8) 


Since  we  need  consider  only  the  space  variation  of  the  potential, 
we  will  henceforth  dispense  with  the  bar  over  the  symbol  for  the 
potential.  In  addition,  for  observation  points  in  the  far  field  of 
the  surface  Equation  (2.8)  simplifies  to: 
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Geometrical  Optics  Result  for  the  Surface  Potential 

To  determine  the  value  of  the  potential  required  by  Equation 
(2.9),  it  is  convenient  to  express  the  monopole  source  at  (0,  h,  0)  in 
a  double  Fourier  series  in  much  the  same  fashion  as  for  instance, 
Brekhovskikh  (1960).  For  y<h,  the  unity  amplitude  velocity  potential 
due  to  the  source  is: 


f0  ■  wr**  i‘v  *  ^  • 


(2.10) 


where 


kly  "  <V  *  k 


lx 


-  k  V/2 

Klz  ; 


and 


R0  =  [x2  +  z2  +  (y-h)2]1/2 


(2.11) 


Since  the  phase  of  the  potential  does  not  change  upon  transmission 
into  the  second  medium,  use  of  the  tangent  plane  approximation  permits 
us  to  write  the  potential  for  the  transmitted  wave,  evaluated  at  the 
point  P  on  the  interface,  as 


00 

>t(x,y0,s)  =  jz  ffi  exp  A[klxx 

_00 


+  klzZ  '  kly(y0-h)])  . 

1  ly 

(2.12) 


where 


2pxc2  cos 

T  “  P2c2cos  +  p^cos  @2 


(2.13) 
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is  the  (plane  wave)  transmission  coefficient  evaluated  at  point  P  on 
the  interface.  Of  course,  we  relate  9 to  9 ^  by  means  of  Snell's 
law;  i.e. , 

kL  Sin  9±  =  k2  Sin  02  .  (2.14) 

To  evaluate  the  potential  at  the  interface,  using  Equation  (2.12), 
apply  the  geometrical  optics  approximation  and  restrict  the  location 
of  the  source  point  such  that  the  quantity  k^R^  remains  large.  In 
this  instance,  the  integral  may  be  evaluated  by  the  method  of 
stationary  phase  in  much  the  same  way  as  that  employed  by  Park  and 
Erteza  (1969).  To  illustrate  the  use  of  this  method,  we  seek  to  find 
the  points  at  which  the  phase  term  in  Equation  (2.12)  takes  on 
stationary  values.  That  is,  defining  the  phase  function  as 


t  =  klj£x  +  k,^z  +  k,„(h-yn)  , 


lz 


ly 


we  seek  values  for  k-  ,  k,  and  k,  at  which 

lx*  ly  lz 


(2.15) 


£i_  =  M_  =  o 

*klx  dklz 


(2.16) 


Using  Equations  (2.11)  and  (2.15),  and  the  condition  (2.16),  the 
points  of  stationary  phase  are  found  as 


k  =  k,  £- 

lx  1 


k  =  r-i  (h-yA) 
ly  RQ  v  JQJ 


lz  1  R0 


(2.17) 
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For  the  geometrical  optics  approximation,  there  is  no  need  to  retain 
terms  in  Equation  (2.19)  above  the  second  order.  It  is  shown  in 
Appendix  A  that,  under  these  conditions,  Equation  (2.19)  may  be 
directly  integrated  with  the  result  that: 
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~  rT  !  ^WVJ0 


V*‘yo’z)  :  tj— ]0  e 

iy 


i[*]"  I  A2*  ,  A2*  ,  A2*  2 

i4"2  0  ‘  [^f^7 lo 

dklx  b\z  lx  lz 


-1/2 


(2.20) 


Using  Equations  (2.15)  and  (2.17),  it  follows  immediately  that 


[♦Jo  =  klR0  • 

r-Sjjfi  W  -  z2> 

9klx2  0  kl  (y0  -  h>2 

»»  <  '  >,) 
L  2  J0  2 

aklz  kl  (y0  •  h) 


(2.21) 


and 


rSLi. 


L  =  - 


R0xz 


Ak  Ak  J0  *  ,  .  .2  ’ 

w  lxu  lz  k1(5r0'h) 

with  which  Equation  (2.20)  may  be  evaluated,  obtaining: 


Vx’yo>z)  ~~  [Tlo  -r 


lkiRo 


(2.22) 


In  an  exactly  analogous  manner ,  the  scattered  potential  at  the  portion 


of  the  interface  adjacent  to  medium  1  can  be  found  to  be  given  as 

(2.23) 


\(x’Vz)  -  1 


,iklR0 


lo 


where 

P2c2cos  t?1  -  p^cos  @2 
S  ~  p2C2COS  ei  +  Dlclcos  S2 


(2.24) 
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is  the  reflection  coefficient  in  which  G 2  is  again  related  to  6 ^  by 
means  of  Snells  law. 


where  the  symbol  <t>t  denotes  the  transmitted  potential  evaluated  very 
close  to  the  point  x,  y^,  z.  Now,  for  surfaces  of  comparatively  gentle 
slope,  the  space  rate  of  change  of  the  transmission  coefficient  may  be 
expected  to  be  small  compared  to  that  of  the  exponential  and 

dk  dk 

Vexp  iCk^  +  k2R2)  — 12  .  (2,26) 

iy 

Moreover,  since  the  maximum  space  rate  of  change  of  the  phase  must 


V*t+(x,y0,z)  -  h. 


occur  along  the  rays, 


V  exp  i(k1RQ  +  k^)  =  i  J  exp  i(k^Rg  +  k^)  ,  (2,27) 

where  kfc  is  the  transmitted  ray  which  is  related  to  the  incident  ray 
by  Snell’s  law.  Thus,  the  gradient  of  the  potential  evaluated  at  the 
surface  becomes: 
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ik 


*t(x’ V‘>  =  ^  I] T  ^  exp(ik1R0)dV.lxdkl2  . 

_00  *"Y 


(2.28) 


Noting  the  similarity  of  Equation  (2.28)  with  Equation  (2.12),  the 
geometrical  optics  solution  to  Equation  (2.28)  follows  immediately  as: 


V*t(x.y0,z)  = 


ikiRo 


i[ktT]o  l 


(2.29) 


Using  a  coordinate  system  centered  at  the  point  P,  ana  new  unit  vectors 

A  A 

a^  and  afc  aligned  with  the  normal  and  tangent  to  the  surface,  it  is 
not  difficult  to  show  that 


<L 

dn 


ikiRo 


ot(x,y0,z)  =  i  k2  [T]q  (-  2 


h-yQ-xa)sH  sin(o)sx-7) 

[1  +  U)g2H2  sin2  (cOgX-7)]1^2 


> 


(2.30) 


where  the  positive  sense  of  the  normal  derivative  is  that  directed 
from  medium  2  into  medium  1. 


Integral  Equations  for  Reflection  and  Transmission 

Substituting  the  expressions  previously  obtained  for  the  potential 
and  its  normal  derivative  into  the  Kirchoff  integral,  we  obtain: 

ikn 

WW  =  4r 


i (klR0  +  k2R2* 
R0R2 


[T], 


[• 


h-H  cos(o)sx-7)-xu)sH  slafogX-r) 
R0U  +  tos2H2  sin2  (o)sx-7)J1/2 


+  lj  ds 


) 


(2.31) 
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where 


2o,c  cos  e. 

[t]  =  — i-2 — —1 - 

0  P2c2cos  e  +  p  c  cos  e : 


cos  0„  = 


c,  2  c„  2 


2  =  [1  .  (_i)  +  (-1)  cos2  0^l/2 


h-Hcos(ci)  x-7)-xm  H  Sin(cc  x-7) 

C08  0,  =  - &  S _ S 

1  -  >■- 2„2  TT  Trrr 


R0[l+a)sH  sin  (oj,  x-7)  ] 


Ro  =  *-x2  +  +  (yQ-h)2] 


2,1/2 


R2  =  [<x2-x)2  +  (Z2-Z)2  +  (y2-y0)2]1/2  .  (2.32) 

In  an  exactly  analogous  manner,  the  reflected  potential  turns  out  to 


ik 

VXl'yl*V=4 - l»\ 


Hcos(a)sx-7)-h+xa)gH  Sin(ajgx-7) 
sin  (a)  x-7)  ]1/2 


where 


(2,33) 


P2c2cos  6»1 
p2C2COS  @1 


p1c1cos  e2 
^1c1cos  e2 


R1  =  [(xrx)  +  (Zl-2)  +  (yry0)2] 


2,1/2 


(2.34) 
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An  expression  which  will  be  useful  in  subsequent  sections  is  that 
expressing  the  ratio  of  the  transmitted  pressure  at  (*2y2Z2^  to  t*ie 
source  pressure  at  unit  reference  distance  and  is  given  by: 


pi 


=  ^  \(X2>W 


(2.35) 


Equation  (2.35)  may  be  readily  evaluated  numerically,  and  the  FORTRAN 
IV  listing  of  a  digital  program  written  to  carry  out  the  computations 
is  included  as  Appendix  B.  The  only  nonzero  contributions  from  the 
integral  arise  from  the  region  of  the  corrugations  and  the  surface 
integral  becomes  a  definite  double  integral  with 


ds  =  [1  +  o)02  H2  sin2  (u^x-?)  ]1/2  dxdz 

J  J 


(2.36) 


Special-Case  Solutions  for  the  Transmission  Integral 

A  closed-form  solution  for  Equation  (2.31)  is  extremely  difficult, 
if  not  impossible,  to  obtain  in  the  general  case.  Even  for  an  extended 
interface,  the  determination  of  the  stationary  phase  points,  which 
enables  one  to  calculate  a  solution  valid  in  the  geometrical  optics 
sense,  leads  to  transcendental  equations  which  can  only  be  solved 
numerically.  Nevertheless,  it  is  worthwhile  to  outline  the  development 
of  such  a  solution  in  general  terms  and  to  evaluate  it  for  special 
cases  in  which  application  of  suitable  parametric  constraints  renders 
the  mathematics  tractable. 

To  compute  the  stationary  phase  points,  we  first  define  the  phase 
function  from  Equation  (2.31)  as. 
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\|r  =  R^Rq  +  ^2^2  *  (2*37) 

then  the  stationary  phase  points,  x^  and  z^,  follow  as  the  real 
solutions  to  the  relation 


fit  _  Si  =  o 

dx  ?z 


(2.38) 


Using  Equation  (2.37)  in  Equation  (2.38),  the  equations  for  the 
stationary  phase  points  may  be  written  as 

xm(1  +  n  r2)  +  yon.Cy2-yom+n  1  *  X2  =  °  '  (2'39) 

om  om 


and 


R2m. 

zm(l-n  F~) 

m  R 

om 


z2  =  o  , 


(2.40) 


where: 


r  2  2  ,  ,  *  2-il/2 

R  =  tXm  +  Zm  +  ^ 

om  m  m  om 


R2m  =  C(x2*Xm>2  +  <VZm)2  +  (y2*yom)2 jl/2 


=  H  cos  (m  x-7)  , 

om  S  m 


om 


yom  =  '  ~  =  “>SH  sin(“SV7)  ’ 


and 


n  =  c2/cx 
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Using  the  stationary  phase  points  obtained  as  the  solutions  to 
the  above  equations,  the  evaluation  of  the  transmitted  potential 
corresponding  to  each  stationary  phase  point  can  be  carried  out  in  the 
limit  of  geometrical  optics  in  a  manner  completely  analogous  to  that 
developed  in  previous  sections  of  this  chapter,  with  the  result  that 


p  cue^^m  k 

=  FI 


m 


ei  Ro¬ 
om  2m 

-  [li-] 


{p2C2+0riBm> 

-1/2 


[*-f]  A 

Sxz  m  bz  m 


dxdz  m 

where,  in  addition  to  the  special  variables  defined  previously,  we 
have: 


(2.41) 


M  =  k  R  +  k  R 
m  1  om  2  2m 


k  =  /[I  +  (y  )2]  1/2  +  hml^oml 

m  \  w  om  R 

om 

2  R„n.2Cl  +  (y'n,)2]  2 

B  =  (1-n2)  -22 - 2S -  +  n2 

m  Lh  -x  y  '  ] 

m  m  om 


h  =  h-y 
m  om 


and 


—  [^4'J  -  5 —  [1  +  h  &  2 

(jq  -v  2  m  R  m  S 

dx  om 


om 


(y  ‘  )' 

om 


(x+h  y  mr 

m  m  om 


om 


] 


+ 


1 


2m 


[l+y2m*S 


2  2 
y  +(y' )'  - 


om 


om 


(x  -x0+yn  y'  Y 
m  2  J  2nr  om7 


] 
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-1  [^2]m  =  f“  (1  *  “ 2>  +  F"  [1  *  -£S“ 

d 2  om  R  2m  *Um 

om  zm 


n2.  HZ  (z  -z  ) 

-£-  [|-f  ]  =  -  _2  (x  +  hv'  )  +  -V51-  (X  -x  +y  y'  )  , 

a)  dxdz  m  3  m  nrom  J  m2  2m7  om 

K  Ka 

om  2m 


y2m  =  y2  '  yom  =  -(d+yom> 

As  an  example  of  the  use  of  Equation  (2.41),  we  proceed  to  solve  the 
comparatively  simple  problem  of  sound  transmission  from  a  mcnopole 
source  to  a  receiver  directly  below  it  when  the  boundary  separating 
the  source  medium  from  the  receiver  medium  is  that  of  an  extended  plane 
(H=0).  This  problem  provides  a  valuable  end-point  check  on  the  theory 
since  a  solution  to  this  problem  using  other  techniques  has  already 
been  reported  by  Weinstein  (1965)  and  Hudimac  (1957).  For  a  receiver 
at  depth,  d,  the  solution  to  Equation  (2.39)  and  (2.40)  is: 


x  =  z  =0 
m  m 


(2.42) 


so  that  Equation  (2.41)  becomes: 

i(k  h+k  d) 

2plC26 

4>t(0,d,0)  -  (p2c2+D1c1)(h+nd) 


(2.43) 


for  the  case  of  a  plane  boundary  and,  of  course,  a  unit  amplitude 
source.  The  expression  for  the  potential  can  be  readily  transformed  to 
an  expression  for  the  pressure  amplitude  at  depth,  d,  in  a  water  medium 
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( 


due  to  the  operation  of  a  one  erg  per  second  point  source  in  the  air 
medium  above  it  and  the  result  is 


I  I  ^P1C1  1 

|pt  <0,d.0)|-  (-J-i)  (h“T 


nd) 


(2.AA) 


This  last  expression  is  precisely  that  obtained  by  Weinstein  (1965)  and 
Hudimac  (1957)  using  other  methods. 

It  is  interesting  to  note  that  the  expression  describing  the  ratio 
transmitted  pressure  to  the  source  pressure  (at  unit  distance),  i.e., 


pt  ^P2C2 

~  (P2c2+pici)  (h+nd> 


(2.45) 


indicates  that  the  transmitted  pressure  level  is  directly  proportional 
to  the  ratio  formed  by  the  characteristic  impedances  of  the  two  media. 
This  proportionality  constant  is  analogous  to  the  transfer  function  of 
a  resistive  voltage  divider  in  elementary  electric  network  theory. 

Another  feature  of  Equation  (2.45)  which  is  noteworthy  is  the 
effect  of  the  velocity  contrast  factor,  n,  as  a  multiplier  on  the 
water  depth  parameter „  The  presence  of  this  term  indicates  that 
changes  in  the  pressure  ratio  brought  about  by  changes  in  receiver 
depth  are  n  times  more  effective  than  changes  in  this  ratio  brought 
about  by  adjusting  the  source  height,,  In  fact,  at  sufficiently  high 
depth-to-height  ratios  for  the  air  to  water  transmission  case,  only 
very  small  changes  in  the  received  pressure  level  will  be  observed  for 
comparatively  large  changes  in  source  height.  It  must  be  kept  in  mind, 
however,  that  Weinstein  (1965)  showed  that  the  velocity  contrast 
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multiplier  is  correct  only  for  the  geometrical  optics  case  for  which 
it  applies;  it  does  not  appear  in  the  limiting  case  of  low  frequencies. 

To  illustrate  how  the  simplified  geometrical  optics  solution  may 
be  evaluated  for  the  more  complex  case  of  a  corrugated  boundary,  a 
digital  program  was  written  to  carry  out  the  analysis  required  to 
locate  the  stationary  phase  points  and  to  subsequently  compute  the 
transmitted  pressure  ratio  by  means  of  the  relation  (2.35)  applied  to 
the  potential  given  by  Equation  (2.41).  The  FORTRAN  IV  listing  for 
this  program  is  included  in  Appendix  C  and  some  examples  of  its  use 
will  be  presented  in  Chapter  IV.  Note  that  when  more  than  one 
stationary  phase  point  is  involved,  the  total  solution  is  the  complex 
sum  of  the  potentials  obtained  for  each  stationary  point. 

To  simplify  the  process  of  finding  th?  stationary  phase  points, 
the  particular  embodiment  of  the  solution  (2.41)  just  described  assumes 
that  z^=0t  so  that  zm=0  from  Equation  (2.40).  The  roots,  x^,  are  then 
computed  using  Equation  (2.39)  under  the  constraint  that  the  receiver 
lies  in  the  x-y  plane  containing  the  source.  Note  that  solution  (2.41) 
is  particularly  useful  when  the  characteristics  **€  sound  transmitted 
through  a  corrugated  interface  of  large  extent  ^re  to  be  found.  In 
such  a  case,  the  numerical  evaluation  of  the  transmission  integral  can 
consume  considerable  computer  time,  whereas  substantially  identical 
results  can  be  obtained  using  the  simpler  solution.  On  the  other  hand, 
numerical  solution  of  the  integral  equation  is  preferred  over  that  of 
the  simpler  solution  represented  by  Equation  (2.41)  when  the  corrugated 
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section  of  the  interface  is  of  only  limited  extent.  In  these  cases, 
one  often  finds  that  only  a  very  few,  perhaps  only  one,  stationary 
phase  point  may  exist.  When  this  happens,  the  approximate  solution 
tends  to  produce  results  that  average  out  the  interference  effects 
arising  from  the  combination  of  various  rays  at  the  receiver  and  thus, 
some  of  the  finer  sound  field  pattern  detail  contained  in  the  integral 
solution  may  be  lo»t.  Examples  of  this  effect  will  be  illustrated  in 
Chapter  IV. 


CHAPTER  III 


EXPERIMENTAL  PROCEDURE 

Description  of  Laboratory  Apparatus 

To  provide  experimental  corroboration  of  the  theory  developed  in 
Chapter  II,  a  complementary  experimental  program  was  carried  out  using 
the  anechoic  tank  facility  at  the  Ordnance  Research  Laboratory  of  The 
Pennsylvania  State  University.  To  simulate  the  monopole  source,  a 
high  quality  enclosed  back  loudspeaker  (tweeter)  was  used  in  an 
unbaffled  configuration.  To  implement  the  finite,  corrugated  boundary, 
several  model  corrugated  surface  sections  were  constructed  of  one- 
eighth-inch  ABS  (acrylonitrile-butadiene-styrene)  sheet  material 
manufactured  by  Marbon  Chemical  Division  of  Borg-Warner  Corp.  Before 
corrugated  models  were  constructed,  however,  the  suitability  of  ABS 
for  acoustically  simulating  an  uneven  water  boundary  was  examined  by 
comparing  data  obtained  on  sound  transmitted  through  a  flat  water 
boundary  to  that  obtained  when  ABS  sheet  was  used  on  the  water  surface. 
Although  these  preliminary  tests  did  show  some  minor  discrepancies, 
the  general  conclusion  was  that  ABS  could  be  used  to  faithfully 
simulate  the  acoustic  properties  of  the  water  surface  when  insonified 
from  the  air  medium,  providing  intimate  contact  is  maintained  between 
the  lower  portion  of  the  sheet  and  the  main  body  of  water. 
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For  the  purposes  of  this  study,  four  model  corrugated  surfaces 
were  constructed  and  photographs  of  these  are  included  as  Figure  2 
through  5,  inclusive.  The  principal  characteristics  of  these  surface 
models  are  summarized  in  Table  1. 

While  in  use,  each  surface  model  is  stabilized  by  lead  ballast 
connected  to  each  of  the  four  corners.  The  ballast  permits  each 
corrugated  surface  to  float  at  a  level  attitude  and  with  the 
corrugations  completely  submerged.  Before  placing  the  surface  model 
on  the  water,  the  corrugated  ABS  sheet  is  treated  with  a  wetting  agent 
to  exclude  air  bubbles  and  ensure  intimate  contact  between  the  water 
and  the  underside  of  the  corrugated  material.  From  the  air  medium, 
this  section  of  the  surface  constitutes  a  stationary,  corrugated 
air/water  boundary.  As  an  additional  precaution  against  reflections, 
all  inside  vertical  surfaces  of  the  model-surface  frames  are  covered 
with  about  an  inch  of  acoustic  damping  material.  Finally,  before  tests 
were  made,  all  other  portions  of  the  tank  surface  were  covered  with 
acoustic  damping  material  to  substantially  attenuate  the  stray 
acoustic  paths  between  the  air  and  water. 

To  facilitate  collection  of  the  laboratory  data,  an  automatic 
hydrophone  positioning  and  recording  system  was  designed,  constructed 
and  installed  in  the  aneohoic  tank.  By  means  of  this  system,  the 
hydrophone  used  to  record  the  underwater  sound  field  could  be 
automatically  swept  over  the  entire  length  of  the  acoustic  tank  and  at 
any  predetermined  depth.  An  electrical  output  derived  from  the  drive 
system  was  used  to  accurately  position  the  pen  of  an  X-Y  recorder  in 
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accordance  with  the  horizontal  position  of  the  hydrophone. 

Since  acoustic  damping  material  was  used  over  all  of  the  tank 
surface  except  that  portion  covered  by  the  model  interface  section, 
direct  observation  of  the  hydrophone  could  not  be  maintained  during 
the  tests.  Consequently,  an  electro-optical  system  was  provided  to 
automatically  shut  down  the  sweep  and  produce  a  spatial  calibration 
marker  to  the  recorder  upon  reaching  a  fixed  position  at  either  end  of 
the  tank. 

Figure  6  shows  the  major  mechanical  portion  of  the  sweep  system. 
The  large  drumlike  assembly  behind  the  tank  was  used  to  adjust  the 
sweep  to  a  given  depth.  The  small  probe  hydrophone  was  mounted  to  the 
center  of  a  horizontal  L  frame  constructed  of  ABS  sheet.  (The  acoustic 
impedance  of  ABS  material  is  very  close  to  that  of  the  water  in  which 
it  is  immersed.)  Also  visible  in  the  photograph  are  the  lights  used 
for  the  optical  position  sensors.  These  sensors,  located  on  the 
traverse,  fore  and  aft  of  the  wheel  assembly,  provide  the  sense  signal 
required  for  automatic  shut  down  of  the  sweep  when  the  leading  light 
aligns  itself  with  a  photocell  mounted  on  the  end  of  the  frame. 

In  Figure  7,  a  block  diagram  of  the  experimental  apparatus,  the 
configuration  depicted  is  one  that  utilizes  narrow-band  random-noise 
(NBRN)  excitation  of  the  loudspeaker  source.  It  was  found  during  the 
course  of  these  tests  that  NBRN  is  the  most  satisfactory  means  for 
obtaining  laboratory  data,  particularly  when  the  complex  patterns 
produced  by  corrugated  interfaces  are  to  be  studied.  The  CW  or 
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warble-tone  (FM)  signals  usually  fail  to  provide  satisfactory  pattern 
definition  because  of  the  existence  of  residual  reflections  in  air 
(such  as  those  that  exist  between  the  speaker  and  the  surface)  and 
underwater  (as  between  the  bottom  and  the  surface).  Supplying  wide¬ 
band  noise  to  the  loudspeaker  (and  filtering  the  received  signal  only) 
is  also  unsatisfactory  because  of  the  power  limitations  of  the  source. 

Calibration 

The  form  of  the  transmission  equations  which  relate  received 
pressure  to  incident  pressure  (or  source  level)  is  particularly 
convenient  for  making  comparisons  with  experimental  data.  Using  these 
equations,  the  theoretical  data  was  computed  as  sound  pressure  in 
decibels  relative  to  the  source  level.  It  could  then  be  readily 
related  to  received  sound  pressure  level  by  adjusting  the  relative 
level  by  the  figure  obtained  for  the  loudspeaker  source  level  (LSL) 
given  by: 

LSL  =  Ggp  +  20  log  E.n  -  29.4  ,  (3.1) 

where 


LSL  =  Loudspeaker  Source  Level 
(dB  re  1  pbar  at  1  ft.) 

G  =  Loudspeaker  (EIA)  Rating 

bp 


(dB  re  0.0002  pbar  and  ln^  at  30  ft.) 


=  Loudspeaker  Input  Level 
(rms  volts) 

In  all  the  experimental  work  to  be  reported  herein,  the  loudspeaker 
input  level  was  maintained  at  +8dBv(rms).  This  figure  represented  the 
maximum  input  voltage  which  could  safely  be  accommodated  by  the 
loudspeaker  without  leading  to  excessive  distortion  or  mechanical 
failure  of  the  voice-coil. 

Calibration  of  the  receiving  system  was  accomplished  simply  by 
recording  the  deflections  produced  at  the  recorder  for  known  input 
signals  and  adjusting  these  figures  to  account  for  the  receiving 
sensitivity  of  the  hydrophone.  Pertinent  calibration  data  obtained  for 
the  loudspeaker  and  hydrophone  are  summarized  in  Table  II. 

Practical  Considerations  in  the  Design  of  the  Experiment 

The  choice  of  the  various  parameters  of  the  corrugated  surface 
models,  as  well  as  the  selection  of  the  excitation  frequencies  used  in 
the  tests,  were  heavily  influenced  by  the  performance  characteristics 
of  the  acustic  source  and  the  anechoic  tank.  In  order  to  develop 
sufficiently  high  source  levels  when  usin^  NBRN,  it  was  necessary  to 
restrict  the  source  center  frequency  to  values  not  in  excess  of  about 
20  kPz.  However,  even  if  this  limitation  did  not  exist,  it  would  be 
difficult  to  employ  source  frequencies  much  higher  than  this  due  to 
the  accuracy  requirements  which  the  use  of  such  short  acoustic 
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Table  II 


SOURCE  AND  RECEIVER  CALIBRATION  DATA 


Frequency 

Gspa 

b 

Sensitivity 

(kHz) 

(dB  re  0.002  pbar 

(dB  re  1  volt 

and  1  ntyj  at  30  ft. ) 

and  1  pibar) 

8 

50.0 

108.2 

9 

47.8 

106.5 

10 

46.5 

105.0 

11 

46.5 

104.5 

12 

46.5 

103.5 

14 

47.5 

102.9 

16 

53.0 

105.5 

18 

51.5 

105.5 

20 

46.0 

105.1 

a 

University  Sphericon  Tweeter 
^Atlant5c  Research  LC10  hydrophone 
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wavelengths  would  impose  on  the  fabrication  of  the  surface  model. 
Instead  of  the  comparatively  simple  heat  forming  process  used  in  the 
construction  of  the  models  used  here,  it  would  be  necessary  to  utilize 
much  more  precise  construction  techniques,  as  well  as  to  provide  much 
more  sophisticated  apparatus  to  measure  and  maintain  the  various 
physical  parameters  used  in  describing  the  experiment. 

The  lower  frequency  limitation  of  the  experimental  procedure  arose 
as  a  result  of  two  considerations.  First,  the  acoustic  absorption  of 
the  tank  liner  was  highest  in  the  frequency  range  near  20  kHz  and 
diminished  drastically  for  frequencies  below  10  kHz.  Consequently, 
any  data  obtained  for  frequencies  below  10  kHz  was  questionable  due  to 
the  reduction  in  received  signal- to-noise  ratio  caused  by  the  increased 
reverberation  in  the  tank.  Secondly,  for  the  longer  wave  lengths,  it 
became  increasingly  more  difficult  to  acoustically  mask  the  unused 
surface  of  the  tank,  and  the  existence  of  stray  acoustic  paths  which 
bypassed  the  region  of  the  corrugated  model  began  to  seriously  degrade 
the  quality  of  the  data.  Thus,  it  was  concluded  that  meaningful  data 
could  not  be  obtained  for  frequencies  less  than  about  10  kHz. 

With  excitation  center  frequencies  limited  to  the  10  to  20  kHz 
ranpe.  the  design  of  the  corrugated  models  was  undertaken  as  a 
compromise  between  requirements  which  tend  to  favor  the  use  of  models 
of  large  physical  size  and  those  which  favor  the  use  of  smaller  ones. 
Whereas  larger  overall  size  permits  a  greater  amount  of  the  incident 
energy  to  be  transmitted  into  the  water  and  also  allows  a  larger  number 
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of  undulations  to  be  used,  larger  models  suffer  from  the  disadvantage 
that  the  receiver  must  be  located  at  a  greater  distance  from  the 
interface  in  order  to  be  in  the  far  field  of  the  extended  radiator  at 
the  surface  (the  corrugated  model).  With  the  depth  of  the  receiving 
hydrophone  limited  to  nine  feet  with  the  sweep  system  installed,  it 
appears  that  the  use  of  surface  models  larger  than  about  one  foot  on  a 
side  could  lead  to  difficulty,  especially  at  the  lower  frequencies. 

In  practice  however,  it  was  found  that  somewhat  larger  corrugated 
models  could  be  used  and  useful  results  obtained  even  in  the  near- field 
of  the  corrugated  model. 


CHAPTER  IV 


DISCUSSION  OF  RESULTS 

Comparison  of  Theory  with  Experiment 

In  this  section,  the  theoretical  data  obtained  by  numerical 
evaluation  of  the  Transmission  Equation  [Eq.  (2.31)]  are  compared  to 
experimental  data  obtained  for  the  case  of  air-to-water  sound 
transmission  using  the  experimental  program  described  in  the  previous 
chapter.  Of  particular  interest  is  the  examination  of  the  accuracy  of 
the  theoretical  predictions  for  sound  transmission  through  a  number  of 
surface  model  corrugations  designed  with  various  waveheights,  slopes, 
corrugation  wavelengths  and  radii  of  curvature. 

The  results  of  the  data  comparison  are  presented  below  for  each 
surface  model  in  turn  beginning  with  model  No.  1.  For  each  model, 
test  results  are  presented  for  various  excitation  frequencies  beginning 
with  the  lowest  one.  This  order  of  presentation  is  then  continued  for 
each  successively  larger  loudspeaker  tilt  angle  used  for  that  particular 
surface  model.  For  loudspeaker  inclinations  different  from  zero,  the 
actual  location  of  the  loudspeaker  is  to  the  left  of  an  imaginary 
vertical  line  passing  through  the  center  of  the  surface  model.  With 
respect  to  the  curves  given,  these  positions  then  correspond  to 
positions  to  the  left  (negative  side)  of  zero. 
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Figure  8  is  a  comparison  of  theory  and  experiment  for  surface 
model  No.  1,  the  largest  of  the  models.  The  data  was  taken  at  a  test 
frequency  of  10  kHz  with  the  axis  of  the  loudspeaker  aligned  with  the 
vertical  line  passing  through  the  center  of  the  surface. 

Figures  9  through  11,  inclusive,  also  represent  data  taken  with 
zero  loudspeaker  tilt  angle,  although  at  higher  frequencies.  For  this 
model  surface,  normal  loudspeaker  incidence  (zero  tilt  angle)  produced 
the  least  accurate  match  between  theory  and  experiment.  Even  so,  the 
data  can  be  seen  to  agree  well  with  respect  to  average  transmitted 
pressure  level,  and  the  differences  which  do  exist  tend  to  be  confined 
to  the  finer  details  of  the  pattern  shapes.  However,  for  the  highest 
frequency  used,  20  kHz,  Figure  11  shows  that  even  the  finer  details  of 
the  pattern  are  in  close  agreement  with  theory. 

The  somewhat  irregular  appearance  of  the  experimental  data  is 
attributed  to  two  factors:  the  fluctuation  in  the  excitation  signal 
which  arises  out  of  the  requirement  to  use  very  narrow  bandwidths  at 
the  loudspeaker  input  to  produce  adequate  radiated  power  levels  with 
NBRN  excitation;  and  the  spatial  irregularities  caused  by  partially 
coherent  reflections.  In  fact,  some  rather  pronounced  excursions  are 
often  associated  with  the  edge  of  the  plots  (the  ends  of  the  acoustic 
tank)  and  particularly  for  cases  in  which  a  null  should  appear  in  the 
data.  This  is  due  to  the  fact  that  equipment  installed  at  the  corners 
of  the  tank  to  support  the  hydrophone  positioning  apparatus  leads  to  a 
higher  incidence  of  reflections  as  the  receiving  hydrophone  approaches 
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the  limits  of  its  sweep.  For  example,  irregularities  in  the  data  which 
appear  to  arise  in  this  way  are  seen  in  Figure  8  in  the  region  to  the 
left  of  the  figure  and  in  Figure  11  at  horizontal  sweep  positions 
between  23  and  43  inches  on  the  right  side  of  the  figure. 

Figures  12  through  16,  inclusive,  are  for  a  loudspeaker  tilt  angle 
of  three  degrees  and  frequencies  in  the  10  kHz  to  20  kHz  range.  In  all 
cases,  agreement  between  average  transmitted  levels  and  fine  pattern 
detail  is  well  within  expected  experimental  limits  for  a  large  region 
containing  the  principal  lobes  of  the  transmitted  sound  field.  Again, 
the  agreement  between  theory  and  experiment  tends  to  improve  with  an 
increase  in  frequency  and  for  the  highest  frequency  20  kHz  (Fig.  16), 
it  is  exceptionally  good  throughout  the  entire  sweep  range.  Figures  17 
through  22,  inclusive,  represent  data  for  a  speaker  tilt  angle  of  six 
degrees.  Agreement  between  theory  and  experiment  is  again  very  good 
for  all  frequencies  in  the  10  kHz  to  20  kHz  range. 

Figures  23,  24  and  25  show  data  taken  for  the  model  1  surface  at 
eight  degree  speaker  tilt  angle.  These  curves,  and  especially  that  of 
Figure  23,  indicate  a  tendency  for  stronger  lobes  to  concentrate  in  the 
region  beneath  the  actual  location  of  the  source  which,  in  this  case, 
is  located  in  the  air  medium  at  about  the  13- inch  horizontal  position. 
This  is  obviously  not  a  general  property  of  the  transmitted  sound  field 
patterns,  however.  At  least  for  small  changes  in  the  loudspeaker 
angle,  some  tendency  for  the  sound  field  pattern  to  retain  its  shape 
but  shift  its  relative  position  in  response  to  changes  in  louspeaker 
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angle  (position)  was  also  noted.  Usually,  the  amount  of  the  pattern 
shift  can  be  estimated  using  Snell's  law.  This  is  to  be  expected 
since  the  portion  of  the  loudspeaker  pattern  which  insonifies  the 
surface  model  is  essentially  omnidirectional.  Hence,  if  the  surface 
section  was  very  large,  one  could  expect  the  basic  underwater  sound 
patterns  to  be  spacially  repetitive  or  periodic  (with  respect  to  a 
vertical  reference  line  through  the  source)  as  the  source  is  moved  in 
a  given  horizontal  direction. 

Figures  26  through  31,  inclusive,  are  examples  of  underwater  sound 
field  data  compiled  for  surface  model  No.  2.  This  corrugated  section 
actually  consists  of  just  a  single  large  wave,  but  has  the  highest 
value  for  the  waveheight  (H)  and  wave  steepness  figure  (2H/X.  )  of  any 
of  the  model  corrugations  used.  The  slope  of  a  tangent  drawn  at  the 
point  of  inflection  of  this  wave  has  an  inclination  of  about  32 
degrees.  Even  so,  it  is  seen  that  satisfactory  agreement  between 
theory  and  experiment  is  again  obtained  for  all  frequencies  an  angles 
of  loudspeaker  inclinations  used.  The  tendency  for  the  pattern  details 
to  agree  more  closely  as  the  frequency  is  increased  is  again  noted  in 
the  data. 

The  satisfactory  agreement  between  theory  and  experiment  for  the 
transmitted  sound  fields  of  models  No.  1  and  No.  2  at  frequencies 
between  10  and  20  kHz  indicates  that  the  theory  is  probably  acceptable 
when  used  for  frequencies  somewhat  below  10  kHz  Unfortunately,  this 
could  not  be  substantiated  with  the  acoustic  tank  used  because  data 


39 


below  10  kHz  was  not  considered  valid  for  the  reasons  previously 
discussed. 

Figure  32  is  an  example  of  a  comparison  between  theory  and 
experiment  for  surface  model  No.  3  at  zero  degree  loudspeaker 
inclination  and  at  a  moderately  high  (16  kHz)  value  of  excitation 
frequency.  The  agreement  between  the  two  curves  is  satisfactory 
although  sufficient  definition  of  the  three  center  lobes  which  should 
be  present  on  the  experimental  curve  is  lacking. 

Generally,  it  was  found  to  be  difficult  to  obtain  experimental 
transmission  patterns  that  were  well  defined  and  symmetrical  for  this 
surface  model  and  for  the  vertical  loudspeaker  position.  Although  all 
the  tests  carried  out  at  12  kHz  and  above  compared  very  well  with 
respect  to  the  transmitted  pressure  levels,  Figure  32  represents  the 
best  pattern  shape  agreement  obtained  during  the  zero  degree  loud¬ 
speaker  tests  It  is  felt  that  mechanical  inaccuracies  in  the 
construction  of  the  model  surface  were  largely  responsible  for  the  lack 
of  pattern  symmetry  observed. 

Figures  33  through  36,  inclusive,  were  obtained  for  a  loudspeaker 
tilt  angle  of  three  degrees.  Comparison  of  the  data  obtained  at  10  kHz 
(Fig.  33)  with  that  obtained  at  12  kHz  (Fig.  34)  illustrates  that 
satisfactory  agreement  between  theory  and  experiment  is  lost  in  the 
vicinity  of  10  kHz  (and  probably  for  frequencies  below  this  value). 

At  these  low  frequencies,  the  experimental  data  is  at  variance  with  the 
theory,  particularly  with  regard  to  the  existence  of  a  broad  null  just 
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to  the  right  of  positions  directly  below  the  center  of  the  surface 
model.  Although  an  experimental  null  in  the  region  of  about  +17  inches 
is  almost  non-existent  at  10  kHz,  a  broad  null  at  about.  30  inches  does 
appear  at  a  frequency  of  12  kHz  and  is  illustrated  by  Figure  34.  Thus, 
it  is  at  frequencies  of  12  kHz  and  above  that  the  experimental  data 
begin  to  display  increasingly  better  agreement  with  theory.  For 
example  at  16  kHz  (Fig.  35)  and  at  20  kHz  (Fig.  36),  both  pattern 
shapes  and  overall  levels  are  in  good  agreement. 

Figures  37  through  41,  inclusive,  represent  model  No.  3  trans¬ 
mitted  pressure  data  at  a  loudspeaker  inclination  of  six  degrees.  As 
in  the  case  of  three  degree  inclination  described  previously,  the  data 
taken  at  10  kHz  disagree  with  respect  to  the  existence  of  a  dip  in 
transmitted  pressure  in  the  region  to  the  right  of  zero  inches. 

However,  in  this  case,  the  discrepancy  is  rather  small  because  the 
pressure  dip  is  not  pronounced.  For  the  data  at  12  kHz  (Fig.  38),  there 
is  a  broad  null  in  the  region  of  +50  to  +70  inches  and,  again,  the 
experimental  data  also  displays  a  null  in  this  region  at  this  higher 
frequency.  At  still  higher  frequencies,  such  as  16  kHz  (Fig.  39), 

18  kHz  (Fig.  40)  and  20  kHz  (Fig.  41),  agreement  between  theory  and 
experiment  once  again  is  very  good. 

The  final  collection  of  experimental  data  to  be  examined  is  that 

which  illustrates  the  transmitted  sound  field  data  for  surface  model 

No.  4.  Referring  to  Table  1,  it  can  be  seen  that  this  surface  model  is 

designed  with  a  value  of  slope  factor  (2HA  )  smaller  than  that  of 

b 

either  models  No.  1  or  No.  2,  but  twice  that  of  model  No.  3.  However, 
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the  wavelength  of  model  No.  4  is  identical  to  that  of  model  No.  3,  so 
that  the  effects  of  increased  slope  and  curvature  on  the  accuracy  of 
the  theory  may  be  examined. 

As  had  been  expected,  the  pronounced  curvature  of  model  No.  4  is 
sufficient  to  produce  a  situation  in  which  the  geometrical  optics 
solution  proves  to  be  deficient  for  predicting  the  experimental  data. 

In  fact,  for  this  surface,  agreement  between  the  data  is  generally  poor 
throughout  the  10  to  20  kHz  frequency  range  and  for  small  speaker  tilt 
angles.  For  example,  Figure  42  illustrates  such  a  comparison  of  data 
at  12  kHz  and  3  degree  tilt  angle.  Although  the  average  predicted 
pressure  levels  in  this  case  tend  to  agree  within  a  few  decibels,  the 
pattern  shapes  display  little  similarity.  The  same  conclusion  is 
reached  for  data  taken  at  higher  frequencies  such  as  that  for  14  kHz 
(Fig.  43),  16  kHz  (Fig.  44),  18  kHz  (Fig.  45)  and  20  kHz  (Fig.  46). 
Furthermore,  from  these  .lata,  it  does  not  appear  that  any  trend  yet 
exists  toward  increased  agreement  between  the  curves  for  higher 
frequencies  although  such  a  trend  is  to  be  expected.  In  fact,  the  data 
at  20  kHz  is  in  particularly  poor  agreement.  It  is  likely  that 
frequencies  much  higher  than  those  available  for  these  tests  would  need 
to  be  used  before  such  a  trend  becomes  discernible. 

Model  No.  4  transmitted  sound  field  data  for  six  degree  loudspeaker 
tilt  angle  is  included  as  Figures  47  through  51.  Although  the  curves 
still  do  not  agree  as  closely  as  they  do  under  similar  circumstances 
using  the  other  model  surfaces,  the  agreement  is  much  improved  over 
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that  obtained  at  three  degrees.  For  the  frequency  range  12  to  18  kHz, 
average  predicted  pressure  levels  are  very  close  to  experimental  data 
and  the  pattern  shape  details  are  generally  acceptable.  However,  the 
data  at  20  kHz  (Fig.  51)  once  again  must  be  considered  one  of  the 
poorest  cases. 

In  addition  to  demonstrating  the  accuracy  of  a  mathematical  model 
for  the  transmission  problem  based  on  a  solution  of  the  geometrical 
optics  type,  the  above  comparisons  between  theory  and  experiment 
provide  useful  insight  into  the  limitations  imposed  by  the  use  of  such 
a  theory.  One  of  the  limiting  conditions  on  the  theory  which  is 
illustrated  by  the  above  work  is  that  imposed  by  the  use  of  acoustic 
wavelengths  which  are  approximately  equal  to  or  greater  than  the 
corrugation  wavelength.  At  such  frequencies,  it  appears  that  detailed 
transmitted  pattern  information  cannot  be  predicted  accurately  by  the 
theory.  Although  the  actual  conditions  under  which  the  theory  begins 
to  break  down  are  rather  vague  due  tc  the  subjective  nature  of  the 
processes  involved  in  matching  pattern  shapes,  it  appears  that  gradual 
loss  of  the  ability  to  predict  pattern  shape  provides  tie  first 
indication  of  theory  deficiency.  Substantial  discrepancies  in 
predicting  average  pressure  levels  tend  to  occur  only  when  tre  acoustic 
wavelengths  are  increased  still  more.  This  is  an  expected  result, 
since  it  is  the  phase  relationships  between  the  transmitted  rays 
reaching  the  receiver  which  largely  establish  the  characteristics  of 
the  transmitted  sound  field  pattern.  Thus,  substantial  variations  in 


the  pattern  shape  details  arise  from  only  relatively  small  discrep¬ 
ancies  in  predicting  these  phase  relationships. 
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As  an  example  of  this  limit  on  the  theory,  consider  the  data 
obtained  for  model  No.  3  at  10  kHz  (Fig.  33  and  Fig.  37).  For  either 
loudspeaker  angle,  noticeable  differences  in  pattern  shape  are  observed 
between  theory  and  experiment  at  10  kHz,  but  these  disappear  when  the 
frequency  is  raised  to  12  kHz  (Fig.  34,  three  degree  speaker  angle)  or 
to  11  kHz  (Fig.  38,  six  degree  speaker  angle).  There  also  appears  to 
be  some  tendency  for  the  data  to  agree  somewhat  more  closely  at  a 
given  frequency  as  the  speaker  inclination  is  increased,  but  inasmuch 
as  the  loudspeaker  pattern  in  the  region  of  the  interface  is 
essentially  omnidirectional,  the  improvement  noted  is  probably  due  to 
the  spatial  relationship  of  the  loudspeaker  to  the  undulations  below 
it  rather  than  on  its  inclination,  per  se. 

It  is  important  to  note  that  these  apparent  deficiencies  in  the 
theory  occur  for  this  model  surface  (Nc.  3)  at  the  frequency  for  which 
the  acoustic  wavelength  in  one  of  the  media  equals  the  corrugation 
wavelength.  Specifically,  Table  1  shows  that  the  wavelength  in  water 
(the  high  velocity  medium)  is  equal  to  the  corrugation  wavelength  at 
10  kHz.  Finally,  we  remark  that  the  tendency  for  the  geometrical 
optics  solution  to  be  deficient  when  the  acoustic  wavelength  in  either 
medium  is  comparable  to  or  greater  tr an  the  corrugation  wavelength  also 
appears  to  be  indicated  by  the  work  of  Asano  (1966)  which  was  cited 
previously.  His  theoretical  data  for  plane  wave  incidence  on  an 
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extended  corrugated  surface  shows  that  the  low-order  coefficients  of 
his  series  solution  for  the  transmitted  potential  tend  to  vary  very 
rapidly  with  frequency  in  the  type  of  situation  described  above,  and 
that  the  range  over  which  these  coefficients  vary  increases  with  the 
velocity  contrast  between  the  two  media.  Evidently,  the  observed 
"instability"  of  the  series  coefficients  relates  to  the  convergence  of 
the  series  solution  and  to  the  tendency  of  the  geometrical  optics 
solution  to  be  deficient  at  lower  frequencies.  The  limitation  on  the 
theory  described  above  may  be  referred  to  e.s  a  lower  bound  on  the 
applicability  of  the  geometrical  optics-type  solution  to  the  trans¬ 
mission  problem.  The  method  can  also  fail  at  frequencies  above  those 
dictated  by  corrugation  wavelength  considerations  alone,  however,  and 
these  attendant  limitations  arise  as  a  consequence  of  other  assumptions 
made  in  the  development  of  the  theory. 

Comparing  the  data  obtained  for  model  No.  4  with  the  results 
obtained  with  the  other  model  surfaces  provides  an  appreciation  cf  the 
effects  of  other  parameters  on  the  applicability  of  the  theory.  For 
this  model  corrugation,  the  predictability  of  pattern  details,  in 
particular,  was  generally  poor  even  up  to  20  kHz,  with  the  exception 
of  data  taken  at  six  degree  loudspeaker  angle  and  frequencies  of  14 
and  16  kHz.  Obviously,  corrugation  wavelength  considerations  alone  do 
not  explain  these  observations  since  the  wavelength  of  model  No.  4  is 
identical  to  that  of  model  No.  3  for  which  accurate  results  were 
obtained  above  10  kHz.  Neither  does  the  influence  of  large  waveheiglt 
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nor  large  slope  account  for  the  difference  since  both  model  surfaces 

No.  I  and  No.  2  have  substantially  higher  waveheights  and  steeper  slope 

figures  than  model  No.  4.  One  feature  which  does  stand  out,  however, 

is  the  fact  that  model  No.  4  has  a  minimum  radius  of  curvature 
2 

=  l/ojg  H)  of  about  one-half  of  that  of  the  other  models.  Thus,  it 
is  likely  that  the  poor  agreement  between  theory  and  experiment  for 
this  model  can  be  attributed  to  the  inapplicability  of  the  tangent 
plane  approximation  for  the  frequencies  used. 

To  test  this  supposition,  we  note  that  Park  and  Erteza  (1969)  give 
Brekhovskikh * s  criterion  for  the  validity  of  the  tangent  plane 
approximation  as: 


4itr  cos  a  »  X  ,  (4. 1) 

c 

where  r  is  the  radius  of  curvature  and  a  is  the  local  angle  of 
c 

incidence.  Evidently,  the  wavelength  (X)  refers  to  that  in  the  high 
velocity  medium. 

The  influence  of  the  local  angle  of  incidence  is  difficult  to 
consider  in  the  present  case,  since  it  takes  on  many  different  values 
along  the  surface  (as  does  r^) .  However,  if  we  assume  that  this  term 
is  equal  to  unity,  we  may  still  obtain  estimates  on  the  limits  of  the 
tangent  plane  approximation  as  follows.  Rewriting  Equation  (4.1)  in 
terms  of  the  excitation  frequency,  f,  we  obtain  the  inequality 
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Using  the  fact  that  the  transmit ted  sound  field  pattern  of  surface  model 
No.  3  can  be  accurately  predicted  by  the  theory  at  frequencies  as  low 
as  10  kHz,  inequality  (4.2)  indicates  that  the  limiting  frequency  for 
the  tangeot  plane  approximation  for  this  surface  is  at  least  of  the 
order  of  15  times  the  characteristic  frequency  factors,  c^Mxr  .  Note 
that,  for  model  No.  3,  the  tangent  plane  approximation  is  valid  for 

frequencies  down  to  the  region  where  X  =  X  ,  but  at  these  frequencies, 

(U  s 

the  tangent  plane  approximation  must  surely  fail  unless  the  wavehelght 
is  oriy  a  small  fraction  of  the  acoustic  wavelength.  Extrapolating 
our  results  for  model  No.  3  to  the  case  of  model  No.  4,  then, 
inequality  (4.2)  indicates  that  the  tangent  plane  approximation  is 
invalid  unless  the  frequency  used  with  model  No.  4  exceeds  20  kHz  and 
this  conclusion  substantiates  our  previous  observations. 

Discussion  of  Experimental  Discrepancies 

During  this  investigation,  supplementary  work  was  carried  out  in 
an  attempt  to  estimate  expected  discrepancies  between  theory  and 
experiment  and  thus  provide  a  reasonable  basis  for  evaluating  the 
accuracy  of  the  solution  techniques.  As  already  indicated,  comparisons 
of  computed  results  with  experiment  indicated  that  generally  good 
agreement  in  both  pattern  shape  and  acoustic  levels  was  obtained  for 
those  cases  in  which  very  low  frequencies  or  very  severe  corrugation 
properties  (small  rc)  were  not  encountered.  It  was  also  discovered, 
however,  that  there  were  cases  for  which  optimum  agreement  between 


47 


theoretical  and  experimental  curves  required  some  minor  adjustment  of 
the  X-axis  spatial-reference  position  between  the  two  curves..  Proper 
adjustment  of  the  zero  reference  then  resulted  in  coincidence  of  the 
two  patterns.  Since  there  was  strong  correlation  between  the  amount 
of  adjustment  required  and  the  date  on  which  various  experiments  were 
carried  out,  it  was  clear  that  these  adjustments  were  required  due  to 
the  presence  of  various  mechanical  errors  in  the  expeximental  setup. 

Some  of  the  possible  errors  that  were  considered  to  be  contrib-' 
uting  co  the  existence  of  this  pattern  dislocation  effect  include,  in 
order  of  importance: 

1.  Errors  in  the  measurement  of  and  uncertainty  in  the 
position,  orientation,  and  alignment  of  the  loudspeaker 
source^,  the  model  surface  section,  and  the  underwater 
receiver. 

2.  Physical  errors  in  the  slope,  size,  and  construction 
of  the  model  corrugations. 

3.  Flexing  and  bending  of  the  model  due  to  stresses 
applied  by  ballast  weight  and  hydrostatic  pressure. 

The  type  of  laboratory  or  theoretical  checks  that  were  made  to 
assess  the  independent  effects  of  position,  orientation,  and  alignment 
errors  included: 

1.  Rotation  of  the  surface  model  by  turning  it  end-for-end 
to  check  for  pattern  variation  when  the  direction  of 
corrugations  are  reversed;  this  evaluates  the  physical 
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symmetry  of  the  corrugated  model  and  provides  a  check 
on  how  accurately  the  surface  level  can  be  established. 

2.  Rotation  of  the  speaker  about  its  axis  to  check  for 
alignment  of  the  mechanical  and  apoustic  axes. 

3.  Movement  of  the  surface  model  to  different  locations 
in  the  tank  to  establish  repeatability  of  established 
test  conditions  and  check  uniformity  of  the  underwater 
acoustic  environment. 

4.  Introduction  of  small  angular  errors  into  speaker  and 
surface  model  orientation  to  observe  the  effect  on  the 
pattern. 

5.  Introduction  of  a  small  positioning  error  into  the 
location  of  the  speaker  relative  to  the  surface  model. 

As  a  result  of  tests  of  the  type  described  above,  the  following 
conclusions  were  reached  regarding  the  relative  importance  of  various 
parameters  on  the  horizontal  positioning  accuracy  of  the  experimental 
underwater  sound  field: 

1.  Rotation  of  the  surface  model  affects  the  structure 

of  the  main  lobes  somewhat  and  can  result  in  uncertainty 
in  the  positioning  of  the  nulls  in  the  data  by  as  much 
as  12  inches.  This  effect  is  complicated  by  its 
dependence  on  speaker  angle,  hydrophone  depth,  frequency, 
and  corrugation  characteristics.  This,  of  course,  is 


true  of  most  of  the  errors  that  were  examined. 
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2.  No  misalignment  of  mechanical  and  acoustic  axes  could 
be  discerned  for  the  loudspeaker. 

3.  Reestablishing  experimental  conditions  demonstrated  that 
repeatability  of  patterns  is  uncertain  wit.  regard  to 
null  positions  by  distances  of  the  order  of  2  inches. 

4.  Errors  in  measuring  the  speaker  and  surface  model 
inclinations  are  very  important.  A  surface  inclination 
uncertainty  of  about  2  degrees  was  observed  to  affect 
the  resulting  underwater  sound  pattern  by  about  10  to 
15  inches  (horizontally)  at  a  9  ft.  depth.  Changes  in 
speaker  tilt  angle  of  about  1  degree,  for  nominal 
inclinations  near  normal,  can  result  in  pattern  shifts  of 
the  order  of  10  to  12  inches  at  a  9  ft.  depth. 

5.  Small  variations  (of  the  order  of  1  inch)  ir.  the 
loudspeaker  height  have  no  appreciable  effect  on 

the  pattern.  However,  inaccuracies  in  the  horizontal 
displacement  of  the  speaker  relative  to  the  center  of 
the  interface  are  significant..  Here,  measurement  errors 
of  as  little  as  1  inch  can  cause  pattern  position  errors 
of  the  order  of  5  to  10  inches. 

There  are  other  sources  of  error  but  these  were  considered  to  b i 
insignificant  and  no  evaluation  of  these  were  therefore  attempted. 
These  include  discrepancies  introduced  by: 
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1.  Approximation  errors  arising  frou  the  use  of  a 
piecewise  planar  approximation  to  the  surface  geometry. 

2.  Neglecting  multiple  scattering  of  sound  at  the 
interface  (In  air  and  in  water). 

3.  Truncation  and  round-off  noise  produced  by  the 
computer  program  during  the  execution  of  the 
arithmetic  operations. 

As  a  result  cf  the  above  work,  it  was  concluded  that  the 
uncertainties  experienced  in  the  substantiation  of  the  computational 
procedure  were  within  reasonable  tolerances,  and  that  experimental 
results  were  generally  in  good  agreement  with  theory  except  in  those 
cases  for  which  the  geometrical  optics  solution  does  not  apply. 

Examples  of  the  Use  of  the  Approximate  Closed-Form  Solution 

In  this  section,  several  examples  are  provided  to  illustrate  the 
use  of  the  approximate  closed-form  solution  derived  previously  and  to 
establish  the  usefulness  of  this  method  as  an  alternative  to  direct 
numerical  integration  of  Equation  (2.31). 

A  comparison  of  the  results  of  the  approximate  solution  with  a 
representative  set  of  experimental  data  is  presented  as  Figure  52. 
Here,  the  experimental  data,  previously  presented  in  Figure  36  (model 
No.  3  at  3  degree  loudspeaker  angle  and  20  kHz),  is  compared  with  the 
results  of  the  approximate  solution.  It  can  be  seen  that  very  good 
agreement  between  theory  and  experiment  has  been  obtained.  Further, 
comparison  of  the  approximate  theoretical  solution  of  Figure  52  with 
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the  theoretical  data  presented  in  Figure  36  shows  that  the  two  methods 
lead  to  very  similar  results. 

For  the  theoretical  data  of  Figure  52,  the  number  of  stationary 
phase  points  found  within  the  interface  boundaries  averaged  about  eight 
and  was  evidently  sufficient  to  produce  adequate  definition  of  the 
shape  of  the  transmitted  sound  field  pattern.  There  are  cases, 
however,  when  tie  location  of  the  source  and  receiver,  as  well  as  the 
interface  size  and  geometry,  lead  to  only  a  very  few  stationary  phase 
points.  For  example,  Figure  53  is  a  comparison  of  theoretical  data 
obtained  for  a  corrugated  surface  with  a  waveheight  equal  to  that  of 
model  No.  4,  but  with  a  larger  wavelength  specification.  Note  that 
these  data  are  computed  for  zero  source  angle  and  a  symmetrical 
interface;  thus,  the  patterns  are  symmetrical  about  the  zero-inch 
reference  point. 

For  eight  foot  source  height  and  nine  foot  receiver  depth 
specifications  used  in  Figure  53,  only  one  stationary  phase  point  was 
found  with  the  receiver  in  the  regions  of  zero  to  47  inches  and  77  to 
110  inches.  For  this  reason,  the  approximate  solution  tends  to  predict 
values  closer  to  the  average  pressure  levels  in  these  regions  and  fails 
to  account  for  the  large  variations  in  the  received  pressure  level 
caused  by  interference  effects.  For  the  region  between  50  inches  and 
75  inches,  three  stationary  phase  points  exist  and  it  can  be  seen  that 
this  leads  to  much  better  prediction  of  the  transmitted  pattern  shape. 
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Figure  54  is  another  comparison  between  the  two  theoretical 
methods,  assuming  the  same  corrugation  properties  previously  used  in 
Figure  53,  but  at  40  kHz  excitation  frequency.  Only  the  region  between 
50  inches  and  75  inches  leads  to  more  than  one  stationary  phase  point 
and  one  again  notes  the  tendency  for  the  approximate  solution  to 
predict  average  levels  rather  than  the  pattern  detail  when  only  one 
stationary  phase  point  exists. 

To  illustrate  the  substantial  improvement  in  the  approximate 
solution  in  cases  for  vMch  a  larger  number  of  stationary  phase  points 
exists,  refer  to  Figure  55.  Here,  the  transmitted  sound  field  is 
again  computed  for  an  interface  with  the  same  corrugation  properties 
used  in  Figures  53  and  54,  but  with  the  source  height  increased  to  20 
feet.  In  this  case,  a  larger  portion  of  the  interface  is  insonified 
by  the  source  and  the  stationary  phase  points  become  more  numerous. 

The  data  is  once  again  symmetrical  about  the  zero- inch  horizontal 
reference  point  and  the  generally  good  agreement  between  the  two 
methods  is  clearly  evident.  Although  the  number  of  stationary  phase 
points  is  about  three,  on  the  average,  only  one  stationary  phase  point 
was  found  for  data  in  the  regions  from  zero  to  10  inches  and  again  for 
the  region  from  50  to  62  inches. 

This  particular  case  (Figure  55)  provides  an  excellent  example  of 
the  utility  of  the  approximate  closed- form  solution  for  predicting  the 
transmitted  sound  field  for  an  extensive  interface.  It  was  possible  to 
obtain  data  for  this  case  with  only  about  two  percent  of  the 
expenditure  of  computer  time  required  by  the  numerical  integration  of 
Equation  (2.31). 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


Statemen  of  the  Problem 

With  the  development  of  increasingly  more  sophisticated  equipment 
to  measure,  record  and  analyze  underwater  acoustic  signals,  it  has 
become  evident  that  the  usefulness  of  underwater  acoustic  sensors  can 
be  enhanced  by  extending  their  range  of  application  to  the  reception 
and  detailed  analysis  of  underwater  acoustic  signals  produced  by 
airborne  sources.  However,  in  order  to  fully  exploit  any  capability  of 
processing  or  otherwise  utilizing  acoustic  energy  transmitted  between 
fluids  in  general,  more  detailed  knowledge  is  required  of  the 
processes  involved  in  the  transmission  of  acoustic  energy  between  fluid 
media  separated  by  arbitrary  boundary  geometries.  To  advance  our 
understanding  of  these  processes,  this  investigation  was  undertaken 
with  the  intention  of  developing  practical  techniques  for  solving  the 
important  basic  problem  of  sound  transmission  from  a  monopole  source 
through  a  finite,  corrugated  boundary  between  fluid  media  and  to 
corroborate  these  techniques  by  means  of  laboratory  data  obtained  for 
the  important  case  of  air-to-water  sound  transmission. 

Procedure  of  the  Investigation 

The  subject  matter  of  this  investigation  is  divided  into  two  major 
sections.  The  first  of  these  deals  with  the  theoretical  aspects  of  the 


study  and  includes: 


1.  The  derivation  of  mathematical  models,  based  on  a 
particular  form  of  the  Kirchoff  radiation  integral, 
which  may  be  used  to  compute  the  properties  of  the 
reflected  and  transmitted  sound  fields  of  a  monopole 
source  after  its  modification  by  a  finite,  uneven, 
fluid- fluid  boundary  composed  of  regular  corrugations 
(sinusoidal  undulations). 

2.  The  development  of  approximate  closed-form  solutions 
to  the  transmission  integral  described  in  (1)  above 
which  are  useful  in  the  study  of  sound  transmission 
through  certain  extended  interfaces  which  are 
insonified  by  high  frequency  sources. 

3.  The  reduction  of  the  equations  involved  in  both  the 
integral  solution  and  the  approximate  closed-form 
solution  (for  sound  transmission)  to  digital  computer 
programs  in  order  to  facilitate  their  numerical 
evaluation. 

4.  Examples  of  the  practical  application  of  the 
mathematical  models  by  computing  the  transmitted 
sound  field  in  numerous  cases  for  which  independent 
evaluation  of  the  quality  of  the  solutions  can  be 


carried  out. 
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The  second  major  section  of  this  study  concerns  the  experimental 
corroboration  of  the  transmission  integral,  and  the  related 
computational  model,  *>v  means  of  experimental  data  obtained  for  the 
case  of  air-to-water  sound  transmission  through  several  finite, 
corrugated  interface  sections* 

Discussion  of  Results 

The  most  significant  results  obtained  during  this  study  were  the 
development  and  experimental  corroboration  of  practical  techniques  for 
predicting  sound  transmission  through  a  finite,  fluid* fluid  interface 
of  sinusoidal  cross  section.  These  tests  of  the  theory  were  carried 
out  by  utilizing  four  laboratory  models  of  corrugated  surface  sections. 
These  model  corrugations  were  designed  to  simulate  different 
corrugation  properties  of  the  interface  and  were  fabricated 
specifically  for  use  in  obtaining  data  on  the  problem  of  sound  trans¬ 
mission  from  air  to  water.  Thus,  by  varying  both  the  orientation  and 
frequency  of  the  loudspeaker  source  used  to  insonify  each  model  surface 
section,  underwater  sound  field  data  for  transmission  through 
corrugated  boundary  surfaces  was  obtained  over  a  comparatively  wide 
range  of  test  parameters  and  for  two  fluids  of  widely  different 
acoustic  impedance.  For  the  first  three  of  the  four  model  surfaces, 
i.e,,  those  designed  with  parameters  within  the  presumed  range  of 
validity  of  the  theory,  very  good  comparison  between  theory  and 
experiment  was  generally  noted  for  all  frequencies  used,  thereby 
corroborating  the  theory. 
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Model  No.  4,  which  was  designed  with  more  severe  corrugation 
properties,  was  employed  to  obtain  information  on  the  limitation  of  the 
theory  when  it  is  applied  to  a  surface  geometry  containing  regions  with 
small  radius  of  curvature.  As  expected,  use  of  this  corrugated  model 
interface  produced  marginal  results  when  the  experimental  data  was 
compared  to  the  predictions  of  the  theory,  particularly  with  regard  to 
the  lack  of  agreement  which  was  displayed  by  the  spacial  details  of  the 
transmitted  sound  field  patterns.  From  these  data,  it  was  shown  that 
the  principal  limitation  of  the  theory  arose  at  low  frequencies  as  a 
result  of  the  inapplicability  of  the  tangent  plane  approximation  used 
in  fitting  the  boundary  conditions  at  the  uneven  interface.  Using  a 
form  of  Brekhovskikh ’ s  criterion  (4.2)  for  the  applicability  of  the 
tangent-plane  approximation,  an  order-of-magnitude  relation  was 
proposed  to  predict  the  lower  frequency  limit  for  the  theory  in  terms 
of  the  velocity  of  propogation  in  the  water  (the  high  velocity  medium) 
and  the  minimum  radius  of  curvature  of  the  corrugated  boundary.  It  was 
shown  that  the  least  lower  bound  on  the  theory  occurred  when  tie 
acoustic  wavelength  in  one  of  the  media  was  equal  to  tbe-  corrugation 
wavelength. 

To  demonstrate  the  use  of  the  approximate  closed -form  solution  to 
the  transmission  integral,  several  exemplary  solutions  obtained  with 
this  technique  were  compared  with  corresponding  ones  obtained  by  direct 
(numerical)  evaluation  of  the  integral.  It  was  shown  that,  even  when 
only  a  very  few  stationary  phase  points  exist,  the  approximate 
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solution  predicted  the  average  transmitted  pressure  level  reasonably 
well.  However »  the  best  results  were  obtained  when  the  interface  was 
large  and  moderately  rough  and  the  source  and  receiver  were  not  too 
close  to  the  interface.  In  such  a  case,  the  stationary  phase  points 
are  sufficiently  numerous  to  permit  a  good  approximation  to  be  made  to 
the  shape  of  the  transmitted  sound  field  pattern  as  well  as  to  the 
average  pressure  level.  Noteworthy  is  the  fact  that  the  approximate 
closed-form  solution  works  well  in  situations  involving  comparatively 
large  interfaces--situations  in  which  numerical  evaluation  of  the 
transmission  integral  requires  a  l?rge  amount  of  computational  effort. 

Conclusions  and  Specific  Contributions  of  This  Study 

As  a  result  of  this  investigation,  theoretical  solutions  have  been 
obtained  for  the  problem  of  transmission  of  sound  from  a  me  nopole 
source  through  a  finite  corrugated  boundary  separating  two  fluid  media. 
These  solutions  have  been  shown  to  be  valid  in  the  limit  of  geometrical 
optics  for  interface  geometries  which  do  not  contain  undulations  with 
a  minimum  radius  of  curvature  of  the  order  of  a  wavelength  of  sound 
(or  less)  in  the  high  velocity  medium. 

The  specific  contributions  of  this  study  can  be  briefly  outlined 
as  follows: 

1.  Integral  expressions*  were  derived  to  predict  the 
transmitted  and  reflected  potentials  due  to  the 
impingement  of  sound  from  a  monopole  source  on  a 
corrugated  interface  separating  two  fluid  media, 
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2.  For  the  case  of  sound  transmission  through  a  corrugated 
interface,  the  integral  expression  for  the  transmitted 
potential  was  reduced  to  an  approximate  closed-form 
solution  valid  in  the  limit  of  geometrical  optics. 

3.  Empirical  data  was  obtained  on  the  transmission  of 
sound  through  several  finite  corrugated  air- water 
boundary  sections  for  various  excitation  frequencies 
and  source  orientations.  It  is  believed  that  this 
is  the  first  time  that  experimental  data  has  been 
reported  on  transmission  of  vibrational  waves  through 
a  rough  boundary  between  fluids. 

4.  Computational  techniques  wen*,  developed  tc  evaluate 
the  transmission  integral  described  in  (1)  above  and 
the  approximate  closed-form  solution  described  in  (2). 

5.  Extensive  comparisons  were  made  between  the  theoretical 
results  obtained  by  evaluating  the  transmission  integral 
and  the  experimental  data  described  in  (4)  above.  These 
comparisons  were  used  to  corroborate  the  theory  and  to 
illustrate  the  conditions  under  which  the  accuracy  of 
the  theoretical  procedures  will  be  maintained. 

6.  Practical  examples  were  provided  to  demonstrate  the 
use  of  the  approximate  closed- form  solution  and 
comparisons  were  made  with  other  data  to  illustrate 
its  suitability  as  an  alternative  tc  numerical 
integration  of  Equation  (2.31). 
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NOT  REPRODUCIBLE 


Figure  3.  Surface  Model  No.  2 
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Nor  RE producible 


Figure  6.  Hydrophone-Sweep  System 


Figure  7.  laboratory  Apparatus  for  Measuring 
Underwater  Sound  Field 


Figure  8.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  1  (Transmitter  Frequency:  10  kHz,  Speak 
Angle:  0  deg.)  K 
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HORIZONTAL  DISTANCE  FROM  CENTER  OF  SURFACE  MODEL  (inches) 

Figure  14.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  1  (Transmitter  Frequency:  14  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  16.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  1  (Transmitter  Frequency:  20  kHz,  Speaker 
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Figure  20.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  1  (Transmitter  Frequency:  16  kHz,  Speaker 
Angle:  6  deg.) 
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Figure  25.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  1  (Transmitter  Frequency:  18  kHz,  Speaker 
Angle:  8  deg.) 
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Figure  26.  Underwater  Sound  Field  Beneath  Surface  Model 
No.  2  (Transmitter  Frequency:  10  kHz,  Speake 
Angle:  3  deg.) 
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Underwater  Sound  Field  Beneath  Surface  Model 
No.  2  (Transmitter  Frequency:  18  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  29.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  2  (Transmitter  Frequency:  10  kHz,  Speaker 
Angle:  6  deg.) 
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Figure  33.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  3  (Transmitter  Frequency:  10  kHz,  Speaker 
Angle:  3  deg.) 
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No.  3  (Transmitter  Frequency:  12  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  35.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  3  (Transmitter  Frequency:  16  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  38.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  3  (Transmitter  Frequency:  11  kHz,  Speaker 
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Figure  40.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  3  (Transmitter  Frequency:  18  kHz,  Speaker 
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No.  3  (Transmitter  Frequency:  20  kHz,  Speaker 
Angle:  6  deg.) 
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Figure  42.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  4  (Transmitter  Frequency:  12  kHz,  Speaker 
Angle:  3  deg.) 
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No.  4  (Transmitter  Frequency:  14  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  44.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  4  (Transmitter  Frequency:  16  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  45.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  4  (Transmitter  Frequency:  18  kHz,  Speaker 
Angle:  3  deg.) 
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Figure  46.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  4  (Transmitter  Frequency:  20  kHz,  Speaker 
Angle:  3  deg.) 
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(  H3AI303d  IV  Joqrf  I  9J  QP)  13A31  3dnSS3yd  QNDOS 


No.  4  (Transmitter  Frequency:  16  kHz,  Speaker 
Angle:  6  deg.) 


NARROW-  BAND  RANDOM- NOISE  EXCITATION 
TRANSMITTER  BANDWIDTH:  50Hz 


(  U3AI303H  iv  I  9J  gp)  13A31  3UnSS3Hd  ONflOS 


Figure  51.  Underwater  Sound  Field  Beneath  Surface  Model 

No.  4  (Transmitter  Frequency:  20  kHz,  Speaker 
Angle:  6  deg.) 
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Figure  53.  Comparison  of  Data  Between  Integral  Equation 
and  Approximate  Closed-Form  Solution 
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Figure  55.  Comparison  of  Data  Between  Integral  Equation 
and  Approximate  Closed-Form  Solution 
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Completing  the  square  within  the  braces  and  rearranging  the  result 
gives: 

00  00 

I  -  /ex p  *1  [(y-B)2  (b-  ~  )]  dy  jexp  [(x-a)  +  ^  (y-B)]2  dx 
•00  .00 

In  the  rightmost  integral,  make  the  substitution: 

u  =  (x-a)  +  |  (y-B) 

and  the  integral  may  be  rewritten  as 
2 

00  iflQ  00 

J  e  2  du  =  2  /[cos  ~  u2  +  i  sin  -  u2]  du 
.00  0 

00 

=  2  (1+i)  /  cos  j  u~  du 
0 

=  (1+i)  (j)1/2  . 

Substituting  into  the  original  equation,  we  have: 

2 

|  [(y-B)2  (b-  j  )]  dy  . 


and  introduce  the  change  of  variables 


i  =  (i+i)  <r>1/2  A*: p 


In  this  equation,  let 


i (b-  r’  ■ r 


u  =  y-B 


Chen 
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i  ‘  (i+i>  <f)l/2 

O 


/ 


oo  9 

i7u 

e  du 


vihich  is 


I  =  2ni  (2«7)*1/2 

or,  in  terms  of  the  original  variables, 

1 — ]  [iv2 
dk 


i  =  2*i  —  .  _ 

2  0..  2  J0  lak,  dk,  J0 


lx 


dk 


.]  -  [<& _ i2 

JA  LnT  \i.  jf 


Iz 


Ix^lz 


{ 


-1/2 
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APPENDIX  B 

COMPUTER  PROGRAM  FOR  NUMERICAL  INTEGRATION 
OF  THE  TRANSMISSION  INTEGRAL 

The  digital  program  Included  in  this  Appendix  was  written  to 
provide  a  means  for  numerically  Integrating  the  transmission  integral 
(2.31).  The  method  employed  is  that  of  carrying  out  the  dcuble 
sunmation  suggested  by  the  integral,  using  small  x  and  y  Intervals 
whose  size  is  set  by  specification  of  the  parameter  DELTA.  The 
successive  contributions  to  the  final  result  are  summed  within  the 
double  DO  loop  using  double  precision  arithmetic.  Within  the  inner  DO 
loop,  the  sign  of  the  quantity  in  the  Snell's  law  radical  is  checked 
and  solutions  representing  negative  values  of  this  quantity  are  cast 
out  since  they  represent  cases  in  which  incident  rays  make  an  angle 
with  the  normal  to  the  surface  which  exceeds  the  critical  angle. 

One  difficulty  which  can  be  foreseen  with  the  use  of  this  version 
of  the  program  is  the  loss  of  accuracy  which  can  occur  when  the  Mray 
path  lengths"  (kjRg+k2R2)  become  very  large.  In  such  cases,  it  is 
well  to  restructure  the  computations  made  to  account  for  relative 
phase  delays  at  the  receiver  in  order  to  avoid  loss  of  precision  in 
summing  the  various  contributions  to  the  magnitude  of  the  total 
transmitted  potential  (or  pressure). 
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NUMERICAL  FVAL'iAT  I  UN  OF  TRANSMISSION  INTFGRAL 


CALCULATION  OF  RATIO  OF  TRANSMITTED  PRESSURE  TO  INCIDENT  PRESSURE 
FOR  THF  SOLUTION  OF  THF  PROBLEM  OF  SOONO  TRANSMITTED  FROM  A 
MONOPOLF  SOURCE  THROUGH  A  CORRUGATED  FIJI  10  RUM  NO  ARY 
J.MACMJtSO  m«l  10/177*9  VERSION  3 


INPUT  UNITS  ARF  TO  BE  SUPPLIED  In  ET*  LRS*ET/SFC*LBS/ET**?tKHZ  f. 
DEGREES.  APPROPRIATE  11*11  Cl INVFRS IONS  ARF  MAOF  BY  THF  PROGRAM. 

LIST  OF  SYMBOLS  FUR  PRINCIPAL  V ARfARLFS 


$YM?hil 


DESCRIPTION  iw  FOII  AT  I  ON 


ANS 

ANSI 

Cl 

C? 

oflta 

FRFO 

Gamma 

h 

NX 

N7 

PHASF 

RHfll 

RHO? 

RO 
R  7 

SRCMT 

T 

WS 

X  I  ,  7. 1 

XSTRT 

V76RO 

X?,V?,7? 

ZSTRT 


WFCFIVFO  PMFSSORF  AMP  RF  SOMRCF  AMP 
?0*LOG| ANSI 

VFLOCITY  OF  PROROGATION-MEDIUM  \ 

VELOCITY  OF  PROP  HO  AT  ION— MFOIMM  7 

X  ANI>  Y  INTFRVAL  SI7F 

FRFOOFNCY  OF  SOMRCF 

CORRUGATION  PHASF  ANGLE 

CORRUGATION  amplitmdf 

NUMBER  OF  GRID  POINTS  IN  X  DIRECTION 

NUMHFR  OF  GRIM  POINTS  IN  7  DIRECTION 

SOMRCF  TO  RFCFIVFR  PHASF  CHAMGF 

DENSITY  OF  MFOIMM  1 

DENSITY  OF  MFOIMm  ? 

DISTANCE  FROM  SOIIRCF  TO  INTFRFACF 
DFSTANCF  FROM  INTFRFACF  TO  RFCFIVFR 
SOURCE  HEIGHT 
PFRIMD  OF  CORRIIGATIMN 
-  TWOPf/T 

VARIABLES  OF  integration 

I  MI T I AL  VALUE  OF  X 
=  H*CnS(WS*Xl-GAMMA) 

RFCFIVFR  COORDINATFS 
INITIAL  VALMF  OF  7. 


R F AL  Kl,«7 
RFAL*«  TMAGl.TMAG? 

NAMFLIST/MATIN/RHIII  ,C1  ,kH0?tC2,  H,  Tf  G  AMM  a  ,  SRCMT ,  DFLT  A ,  XSTR  T  ,  l  STRT , 
lFREOfNX,M7/RCVRPT/X7tY?f 7?, X?lNfY?IN, ZPIN/lim /ANS, AMS  1 , TMAG1 »TMAG2 
?,FPHSF 

l  RFAD( 5,2 )  RHOl f Cl ,WH02,C?,H, T, GAMMA, SRCHT,DFLTA, XSTRT, 7 STRT, FRFO, 
1NX,NZ 

?  FORMAT! 1 2F5.1 ,2IS) 

WR  1  TF  ( f»*OA T !  N  ) 

TW0PI=7, *3. 141 5926 
FRFOsFRFO«ir»no, 

WSa  TWOPI/T 

GAMMA=  GAMMA*  {  TWflPI  /  ) 

C R A  T 1 0  =  C2/C1 
^  RFAD(5,4)X?,Y7t7? 
i*  FORMAT  (  3Fi  n.3  ) 

X2 I N=  X  7* 1 7, 

Y?1N  =  Y7*  )?. 

7?\N=7?+\ 7. 


WRI  TF(6f  Rf.VRPT ) 

*»  IF(Y2.F0.n,  )  STOP 
K\m  fTWnPI*FRrO)/C1 
K?»  (TWflP|*FRFO)/C? 

4  TMAG1=0. 

TMAG2*0. 

7  nr»  «  ?«i,nx 

XI*  XSTRT  *  I*nFLTA-nFl.TA/?. 

YZERO  *  H*f,MS(WS*XI-GA*MA) 

nYZRn*  ws*h*si N( ws*x i -gamma  j 
FCTRl*  {  «;Rr.MT-Y7FRM-Xl*nYZRfl) 

FCTRIP  *  E0TR1**? 

nn  «  .1*1  *N7 

71*  ZSTRT  +  .l#OFLTA-DFL  TA  /  ?. 

R02*  (  X I **?*7 !**?♦(  YZERn-SRf.HT)**?) 

R2  *  SORT  (  ( X  ?-X I  )**?  +  (  7  2-7. 1  )*#?+(  Y 2- YZERO )  **?  J 
RO* SORT(  R rt?) 

FCTR2*  »n?*n .+nY7»o**;>) 

RAD1CL*  (1  .-r.RATin**2)*(  PCTR?/FCTR1P)+  CRATIO*#? 

IFIRADICL.LT.O.  )GO  TO  H 

FCTR3*  RHn?*C?+  RHtM*C l* SORT (RAO? CL  J 

FCTR4*  ( ^CTRl /( RO^SQRT ( l •♦OYZRO**2 ) ) )♦  1.0 

OELMG*  (  |FCTR4*(DFLTA**?  ) ) / ( R0*R?*FCTR3 )  )*SQRT(  1  .+OYZRO**?) 

PHASE  *  K | $MO+  K2*K2 

phase*  AMrtn<, phase*  two  n 

TMAGl*  TMAG1+  OELM(;#Cf3S  (  PHASE  ) 

TMAG2*  TMAG?  +  DELmG*S1N( PHASE) 

R  CONTINUE 

EPHSE*  TMAGP/TMAG1 
EPHSE*  AT AN( EPHSE ) 

EPHSE*  AMRO ( EPHSE  *  Tw0 P 1 ) 

EPHSE*EPHSE*(3*0./rw0P) ) 

ANS*  RHn?«EREO*nSORT( TMAG1**2+TMAG?**?J 
ANSI*  (  Al.nGim  ANS))»20. 

WRITE! 6*  OUT ) 

GO  TO  3 
FND  ' 
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APPENDIX  C 

COMPUTER  PROGRAM  FOR  EVALUATION  OF 
APPROXIMATE  CLOSED-FORM  SOLUTION  OF 
TRANSMISSION  INTEGRAL 

The  digital  program  included  in  this  Appendix  was  written  to 

facilitate  evaluation  of  the  roots  of  the  stationary  phase  relations, 

Equation  (2.39)  and  (2.40),  for  the  specific  case  z^O  and  thereby 

to  compute  the  transmitted  pressure.  This  particular  program  obtains 

the  stationary  phase  points  in  two  stages  or  steps.  First,  an 

estimate  of  the  locations  of  the  roots  is  found  by  checking  the  sign 

of  the  result  for  Equation  (2.39)  as  x  is  varied  from  a  starting 

m 

value,  XSTRT ,  to  its  final  value,  XSTOP,  at  intervals  of  DELTAX. 

Possible  real  roots  are  indicated,  and  their  values  stored  in  vector 

ROOTS(I),  when  particular  successive  values  of  x  cause  the  left-hand 

m 

side  of  Equation  (2.39)  to  vanish  or  to  alternate  in  sign.  After  all 
the  possible  approximate  roots  are  catalogued  in  this  way,  these  roots 
are  more  closely  approximated  by  the  use  of  subroutine  RTWI  as  the 
second  step  in  their  evaluation.  After  locating  all  the  points  of 
stationary  phase,  the  program  evaluates  the  total  potential  at  the 
observation  point  by  summing  the  individual  contributions 

using  Equation  (2.41).  Finally,  the  tot^l  transmitted  pressure  is 
calculated  in  decibels  and  presented  as  output. 
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EVALUATION  OF  POINT-SOURCE  CORRUGATED  INTERFACE  SOLUTION  BY 
GEOMETRICAL  OPTICS  APPROXIMATION  OF  TRANSMISSION  INTEGRAL 
SURFACE  FOUATION  IS  YO«H*CQS(WS*X-GAMMA) 

SOURCE  POSITIONED  DIRECTLY  OVER  ORIGIN 
RECEIVER  IN  X-Y  PLANE 

J.MACALUSO  ORL  11/20/69  VERSION  3  MOD  1  11/22/69 


LIST  OF  SYMBOLS  FOR  PRINCIPAL  VARIABLES 


SYMBOL 


DESCRIPTION  OR  EQUATION 


Cl 

C2 

0 

OFLTAX 

F 


VELOCITY  OF  MEDIUM  1 
VELOCITY  OF  MEDIUM  2 
RECEIVER  DEPTH 

INCREMENT  OF  X  FOR  ROOT  SEARCH 
FRFQUFNCY 


GAM 

GAMMA 

H 

ITER 

Kl 

K2 

NN 

PHASE 

PHI  TOT 

PHIRE 

PHI  |M 

RHOl 

RH02 

ROM 

ROOT!  I  ) 
R2M 
SRCHT 
T 

WS 

XMP 

XM 

XSTRT 

XSTOP 

Y7ERO 


CORRUGATION  PHASE  ANGLE! DEG) 
CORRUGATION  PHASE  ANGLE (RAO I 
HEIGHT  OF  CORRUGATION 
NUMBER  OF  ITERATIONS  IN  RTWI 
WAVENUMBER  IN  MEDIUM  1 
WAVENUMBER  IN  MEDIUM  2 
C2/C1 

SOURCE  TO  REC  PHSE  ANGLE  FOR  XM 
TOTAL  (REL)  PRESSURE  AT  RECEIVER 
INTERMEDIATE  VALUE  OF  POTENT  I AL ( RE  PT ) 
INTERMEDIATE  VALUE  OF  POTENT  I AL < IM  PT) 
DENSITY  OF  MEDIUM  1 
DENSITY  OF  MEDIUM  2 
DISTANCE  FROM  SOURCE  TO  STAT  PHSE  PT 
I  1 TH  APPROXIMATE  ROOT 

DISTANCE  FROM  RECEIVER  TO  STAT  PHSE  PT 

SOURCE  HEIGHT 

PERIOD  OF  CORRUGATION 

?*PI/T 

TRIAL  VALUE  OF  X  IN  ROOT  SEARCH 
MOST  ACCURATE  VALUF  FOR  I • TH  ROOT 
STARTING  VALUE  OF  X  FOR  ROOT  SEARCH 
FINAL  VALUE  OF  X  FOR  ROOT  SEARCH 
VALUF  OF  Y  AT  COORDINATE  XM 


DIMENSION  ROOTS  I  50  ) 

REAL  K l ,K2, NN, NN2 

COMMON  NN fWSfH, GAMMA, SRCHT tD,X? 

EXTERNAL  FCT 

NAMELlST/DATIN/RH0UClfRHn2»C2,TtH, SRCHT, D,F,  GAM  ,  i  TER  , XSTRT, 
1XST0P,DELTAX,X2 

1  READ(5,2)RH01,C1*RH02,C2,T,h, SRCHT, F, GAM  , I TFR 

2  FORMAT! 9FB, 1,215) 

3  IF(RH01,E0„0. )STnP 
PI*3. 1415926 
GAMMA»GAM*(PI/lBO. ) 

301  READ ( 5, 300 )  0, XSTRT , XSTOP, DFLTAX ,  X2 
300  FORMAT ! 5E10,3) 

WR ITF ! 6, OAT  I N ) 

I F( D,EQ,0* )  GO  TO  l 
TWOPI*  2*#PT 
WSoTWOPl /T 
B«(WS**2)*H 

Kl*  (TW0PI*F*1000. >/Cl 
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K2*  t^WOP^FMOOO.l/CZ 

NN.C2/C1 

NN2*NN'*2 

PHIRE.a* 

PHI IM*0» 

IOX*0 

NX*(XST0P-XSTRT)/0FLTAX  M, 

XMP.XSTRT-OFtTAX 

00  100  1*1, NX 
XMP-XMP+OFLTAX 
XMPPl-XMP+OELTAX 
V 1*FCT ( XHP ) 

Y2*FCT ( XMPPl ) 

1FUXMP-YU*(XHPP1-Y2)  1101,101,100 
101  IDXaIDX+1 

ROOTS!  !0X)*XMP*0EL  TAX/?. 

IF!  I0X.EO, l  )Gll  TO  100 

1F(  ROOTS !  I  OX ) -ROOT  St  IDX-1)#GT*1  ,5*0ELTAX )G0  TO  100 
inx*mx-i 

ROOTS! IOX)*ROOTS< I  OX ) ♦ORLTAX/2. 

100  CONTINUE 

00  105  KK*1,I0X 
105  WRITE!6,10M  KK,R00TS(KK ) 

I0f>  FORMAT!  1H  ,  1 3, •  APPROX  1  MATE  ROOT  FOUND  AT  XM*»,F10.3) 

no  r  j*i,  mx 
XMP-ROOTS!  .1) 

CALL  RTWl(XMtvAL,FCT,XMP, l.E-06, ITER, !ER ) 

I ER*  1 ER+ l 

GO  TO (501, 502, 503), l ER 
50?  WRITE! 6»5l?)XM,XMP 

51?  FORMAT! 1H0, »N0  CONVERGENCE,  XM*» , E 10. 3, ' XMP* ' , E 10.3 ) 

GO  TO  410 
503  WRITE(6,513) 

513  FORMAT ( 1  HO, ' OENOMl NAT0R*0# , I #E*  FNCN  OER IVAT IVE* 1 , CANNOT  CONTINUE 
1G0  TO  NEXT  RECEIVER  PO*NTM 
GO  TO  301 

501  WRITE(6»51 1 )  VAL,XM,XME 

511  FORMAT ( 1H  , 1 VAL  OF  X-F(X)  AFTER  CONVERGENCE* ',610.3,' XM* ' ,610.3, 

1 »XMP**,E10.3) 

HT*SRCHT-H*COS(WS*XM  -GAMMA) 

411  XM2«XM**2 

nn»04H*cns(ws*xM-GAMMA i 
R0M2*  XM2+HT**2 
R2M2*( X2-XM)**240D**2 
ROM*SORT( R0M2 1 
R2M»SORT(R?M2) 

YZERO*H*CnS(W.S*XM-GAMMA) 

YP»WS*H*S1N<WS*XM-GAMMA) 

YP2»YP**? 

RAOICL*  (  !  (1,-NN2  )*ROM2*n.+YP2) )/(  (  HT-XM*YP  )**2 )  )  *NN2 

IF(RADICL.GE,0.1  GO  TO  11 
410  PHIT0M*0, 

GO  TO  12 

11  RADICL*  SORT! RAOICL  1 

FCTR1*  RHO?*C2+RH01 *C1*RA0ICL 

FCTR2*  (1.4(  HT-XM*YP)/(R0M4S0RT( U4YP2) ))*SORT( 1.+YP2) 

FCTR3* ( NN/ROM )* ( 1 ,+HT*YZERO* ( WS**2 1+YP2-! ( XM+H1 *YP 1**2 ) /R0M2 ) 
FCTR4* ( 1 .-0D*( WS**2 14YZFR04YP2-! ( XM-X2-nD*YP 1**2 1 /R2M2 ) /R2M 
02BDX*  FCTR34FCTR4 
02B0Z  *  NN/ROM  ♦  1./R2M 


FCTR5  ■  02Rnx*02BDZ 
FCTR5*  ABSIFCTR5) 

FCTR5*  SORT ( FCTR5) 

PH1TOM  -  (RH01*C2*FCTR2)/(FCTR1*FCTR5*R0M*R2M) 

9  PHASE  ■  Kl*Rf1M4-K2*R2M 
PHASE  >  AMnO( PHASE *TWOPI) 

PHIRE-  PH! TOM*COS( PHASE)  +  PHJRE 
PHI  IH*  PHlTOM*SIN(  PHASF  )  ♦  PHIJM 
1?  PHITOT*  SORT<PH|RE**2+PHI IM**2) 

!F(PH1T0T.NE.0.)G0  TO  10 
PHIT0T*-9»  99E+09 
GO  TO  401 

10  PH1T0T  «  ?0.*AL0G10( (RH02*PHIT0T)/RH01) 

401  WRITE<6,4)  JtPHITOT 

6  FORMAT! 1H0, *REL  PRESS  LVL  WITH*, 15, ‘ORDER  CORRECTION  TERM  IS* 
1E10*3, 'OR* /// ) 

IF(H.E0,0.)  GO  TO  301 
A  CONTINUE 
GO  TO  301 
END 

REAL  FUNCTION  FCT1X) 

REAL  K1,K2«NN,NN2 

COMMON  NN,WS*H, GAMMA, SRCHT,0,X2 

HPRME-  SRCHT-H*COS(WS*X-GAMMA) 

OPRME— D-H*COS!WS*X-GAMMA.) 

ROMP*SORT ( X**2+HPRMF** 2 ) 

R2MP«S0RT( (X2-X)**2+OPRME**2) 

YPRME«WS*H*SIN(WS*X-GAMMA) 

FCTM«-(NN*HPRME*YPRME*R2MP)/(NN*R2MP+ROMP) 

FCTR2* ( DPRMF*YPRME*ROMP )  / ( NN4R2MP+R0MP ) 
FCTR3«!X2*ROMP)/(NN*R2MP+ROMP) 

FCT*FCTR 1-FCTR2  +FCTR3 

RETURN 

END 
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